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1 Introduction 



Statistical physics is fundamental physical theory, which deals with equilibrium (or even non- 
equilibrium) properties of a large number of particles using the well established concept based 
either on classical or quantum mechanics. With respect to this, the term statistical mechanics 
is often used as a synonym to statistical physics that covers probabilistic (statistical) approach 
to classical or quantum mechanics concerning with many-particle systems. The most important 
benefit resulting from this theory consists in that it relates microscopic properties of individ- 
ual particles to observable macroscopic (bulk) properties of matter. Even although relations 
between some macroscopic properties and fundamental properties of individual particles are oc- 
casionally elementary (for instance the total mass is simply a sum over particle masses), many 
material properties cannot be simply elucidated from the fundamental properties of constituent 
particles, i.e., from the microscopic point of view. In particular, statistical mechanics enables 
to explain observable macroscopic features of real materials solely by imposing forces between 
the constituent particles. For this purpose, one necessarily needs just some plausible assump- 
tion about internal forces between constituent particles (inter-particle interactions) in order to 
make relevant theoretical predictions for observable properties of a given macroscopic system. 
This assumption, which is built on some realistic microscopic idea of how individual particles 
interact among themselves, constitutes a framework for some simple theoretical idealization to 
be referred to as a statistical model. 

Of course, each statistical model serves only as an approximative description of physical real- 
ity aimed at describing observable macroscopic properties preferably quantitatively or leastwise 
qualitatively. However, it is very difficult and often incredible task to define a realistic model, 
which is on the one hand mathematically tractable and on the other hand provides a compre- 
hensive description of all observable macroscopic properties. The most formidable difficulties 
are usually encountered when attempting to formulate and to solve the relevant model mathe- 
matically. There are just few valuable exceptions. The most common example surely represents 
an exactly solvable model of an ideal gas (no matter whether consisting of classical particles or 
fermions or bosons) in which the constituent particles do not interact among themselves until 
they undergo perfectly elastic collisions. If the inter-particle interactions are taken into account 
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(suppose for instance the real gas instead of the ideal gas), however, realistic models are highly 
appreciated if they are still exactly solvable, but this is usually not the case. 

If the constituent particles of some interacting many-particle system are situated on discrete 
sites of a crystal lattice and only short-ranged inter-particle interactions need to be considered, 
then a substantial simplification in the mathematical treatment of relevant model (s) is usually 
achieved. Under these simplifying constraints, one concerns with so-called lattice- statistical 
models that are generally more amenable to an exact analytical treatment even though so- 
phisticated mathematical methods must be still employed for obtaining exact solutions of even 
relatively simple-minded models. Hence, it follows that exactly solved models are usually con- 
sidered as an inspiring research field to emerge in the statistical mechanics, which regrettably 
requires a considerable knowledge of sophisticated mathematics and are therefore beyond the 
scope of standard courses on the statistical physics. Apart from this drawback, the topic exactly 
solvable models in statistical mechanics surely represent an exciting research field in its own 
right as convincingly evidenced by a rather rich and instantly growing list of excellent books 
devoted to this intriguing subject matter [THT3]. 

The main goal of this book is to make a brief introduction into the method of algebraic 
mapping transformations, i.e. an exact mathematical technique, which enables after relatively 
modest calculation a rigorous analytical treatment of diverse more complex lattice-statistical 
models by establishing a precise mapping correspondence with simpler exactly solved lattice- 
statistical models. 
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2 Ising and Heisenberg models 

In this section, let make few comments on two generic lattice-statistical models, which are of 
particular research interest because of their usefulness and flexibility in representing diverse 
real-world systems of very different nature. 



2.1 Ising model 

The Ising model perhaps represent the most versatile model of statistical mechanics at all, 
which is simultaneously fully mathematically tractable on one- and two-dimensional (ID and 
2D) lattices. This simple-minded lattice-statistical model has been proposed by Lenz in 1920 
[H] and five years later has been exactly solved by Ising for the particular case of the linear 
chain [T5]. Note furthermore that there are several excellent review articles on the historical 
developments of the Ising model to which the interested reader is referred to for further details 

[IMS]. 

The spin- 1/2 Ising model can be defined on a crystal latticd_| through the Hamiltonian 

N 

H = -J (2.1) 

(ij) »=1 

where cr, = ±1/2 is two- valued Ising spin variable situated at the zth site of a crystal lattice, the 
former summation takes into account a configurational energy associated with the interaction 
between the nearest-neighbour spins and the latter summation accounts for the Zeeman's energy 
H = giisB of magnetic moments in an external magnetic field B (g stands for Lande g- factor 
and fiB is Bohr magneton). At first sight, the Ising model might seem to be the greatly oversim- 
plified model as it first takes into consideration only extremely short-ranged interaction^ and 
second, it also neglects all quantum effects in that it disregards a quantum-mechanical nature 
of spin by considering spin as a classical two-valued variable. While the former restriction is 
rather well satisfied in a variety of insulating magnetic materials, the latter restriction turns out 
to be much more profound as far as the theoretical modeling of insulating magnetic materials 
is concerned. From this point of view, the Ising model offers merely a semi-classical description 



1 Note that the Ising model can be defined for systems without translational invariance as well. 
2 The first summation is usually restricted just to the pairs of nearest-neighbour spins. 
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of interacting many-particle systems, because a set of discrete spin values is the only quantum 
feature of this lattice-statistical model. 

Before proceeding to a survey of the exactly solved Ising models, let us briefly comment on 
their possible experimental realizations. The Ising model was for many years merely regarded 
as the purely academic model without any correspondence to a specific real- world system, since 
the first insulating magnetic materials that would satisfy its specific requirements have been 
discovered almost a half century after its invention. At present, there are two wide families of 
insulating magnetic materials whose magnetic behaviour is generally in accord with theoreti- 
cal predictions of the Ising model. The first class of the Ising-like magnetic materials involve 
rare-earth compounds such as Dy(C2H 5 S04)3.9H 2 0, Dy 3 Al 5 0i2, DyPO/t, LiHoF 4 , LiTbF 4 and 
so on (see Ref. [20] and references cited therein). In this class of materials, the only carri- 
ers of magnetic moments are the rare-earth elements (Dy, Ho, Tb, etc.) that interact among 
themselves almost exclusively through dipolar forces, while other non-dipolar interactions are 
in general negligible. Under certain conditions, the dipole-dipole interaction between the more 
distant magnetic moments can be ignored and it is sufficient to consider merely the interaction 
between the nearest-neighbour magnetic moments (Ising-like criterion). It should be never- 
theless pointed out that the magnetic dipole-dipole interaction is long-range interaction in its 
character (it decays as an inverse third power of the distance) and hence, the interactions with 
more distant rare-earth ions can occasionally be at an origin of more complex behaviour. Thus, 
the most suitable experimental realizations of the Ising-like models provide rare-earth com- 
pounds in which the interactions between more distant magnetic moments almost completely 
cancel out and the interaction between the nearest-neighbour magnetic moments makes the 
most profound contribution to the overall magnetic behaviour. 

The second class of the Ising-like magnetic materials represent the insulating magnetic ma- 
terials from the family of molecular-based compounds [2T]. In this wide family of magnetic 
compounds, the most important interaction between magnetic centers (mostly transition-metal 
ions) is a pairwise spin-spin interaction mediated via intervening non-magnetic atom(s) through 
the indirect superexchange mechanism [22H21]- Even though there does not exist a general 
theory, which would admit a straightforward calculation of the interaction parameter J orig- 
inating from the superexchange mechanism and this parameter is usually obtained only as a 
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self-adjustable parameter from a comparison with relevant experimental data, a strength of the 
spin-spin exchange interaction decays very rapidly with a distance - at least as an inverse tenth 
power of distance [25]. From this point of view, the molecular-based compounds involving mag- 
netic metal centers satisfy much better the necessary criterion for the Ising-like materials, which 
demands a predominant nearest-neighbour interaction and negligible interactions between the 
more distant spins. The most crucial limitation, which prevents the most of molecular-based 
compounds to be good examples of the Ising-like materials, thus lies in a demand of having the 
extremely anisotropic exchange interaction while the superexchange mechanism gives rise to 
the isotropic exchange interaction [2T] . However, the anisotropic exchange interaction need not 
arise from the interaction mechanism alone, but it may have a close connection with another 
sources of the magnetic anisotropy such as spin-orbit coupling, crystal- field effect, dipolar inter- 
actions, etc [21HI21]. Under these circumstances, the theoretical description based on the Ising 
model is justified even if it still represents a certain oversimplification of the physical reality. It 
is noteworthy that the overall agreement between theoretical predictions derived from the Ising 
model and the relevant experimental data is generally found to be very satisfactory mainly for 
several cobalt-based compounds such as Co (pyridine) 2 C1 2 , A 2 CoF 4 (A = K, Rb), Cs 3 CoX 5 (X 
= F, CI, Br) [20l[2TJ. 

Finally, it should be also remarked that all magnetic compounds from both the aforemen- 
tioned families of the Ising-like materials are three-dimensional crystals in reality, however, some 
of them can effectively possess the low- dimensional magnetic structure on behalf of the lack of 
an appreciable magnetic interaction in one or more spatial directions. As a matter of fact, the 
magnetic and crystallographic lattices may significantly differ one from each other especially 
if the carriers of magnetic moment are well separated along some spatial direction(s). Conse- 
quently, the magnetic lattice then becomes low- dimensional due to the short-range character of 
magnetic forces and it is therefore of fundamental importance to investigate low-dimensional 
spin models as well. The reliability of exactly solved low- dimensional Ising models in repre- 
senting real-world insulating magnetic materials has been checked with an appreciable success 
even if some healthy skepticism is always appropriate if one is seeking true understanding of 
real magnetic materials [2"U] . 

In conclusion, let us also briefly mention other possible (non-magnetic) applications of the 
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Ising model and its different generalizations like the Blume-Capel model [26} 12?]. the Blume- 
Emery-Griffiths model [2S] and others, which provide a deeper understanding of cooperative 
behaviour inherent to interacting many-particle systems from seemingly diverse research areas. 
Even though the Ising model has been originally invented for describing cooperative nature of 
spontaneous long-range order to emerge in the magnetic materials, throughout the years it has 
proved its usefulness by investigating the order-disorder phenomena in metal alloys [29-3 1J, 
the vapour-liquid coexistence curves at liquid-gas transition |29 ]l30ll32] . the phase separation in 
liquid mixtures [28l|32l[33], the saturation curve of hemoglobine [SIHSHHS], the initial reaction 
rate of allosteric enzymes jHSlEH], the melting curve of helix-coil transition in DNA [3T1 - HD] . the 
role of socio-economic interactions in determining business confidence indicators [HH13], urban 
segregation the language change [13JII7JII8] and many other topics. 

2.2 Survey of exactly solved Ising models 

A great deal of research interest aimed at searching various exactly solvable Ising models has 
resulted in a rather extensive list of up to date available literature concerned with the exactly 
solved Ising models. It is therefore beyond the scope of this work to review all of them and 
selected examples should mainly serve only for illustration and are chosen so as to reflect 
author's previous and current research interests. It should be also emphasized, moreover, that 
the exact solution of any non-planar Ising model is essentially NP-complete problem as recently 
pointed out by Istrail [49J and Cipra [50J, which implies that any 3D Ising model is possibly 
analytically intractable problem despite its conceptual simplicity. Hence, the subsequent list 
of exactly solved Ising models essentially contains only the lattice-statistical models defined on 
ID and 2D lattices. 

It has been already mentioned previously that the spin- 1/2 Ising model on the linear chain 
has exactly been solved by Ising in 1925 with the help of combinatorial approach [T5]. Af- 
terwards, the Ising's exact results have been re-derived using a variety of other mathematical 
techniques such as the transfer- matrix method due to Kramers and Wannier [51]. Among the 
most interesting rigorously solved ID Ising models one could mention the spin-S* Ising linear 
chain [52TI51] . the spin- 1/2 Ising chain accounting for the interactions between the more distant 
spins [55T457] . the spin-1/2 Ising model defined for two, three or four coupled chains [58-62J, 
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the spin-1/2 Ising model with alternating bonds [63] and the mixed-spin Ising chains [641469] 
and ladders [THIEE] supplemented by the zero-field splitting parameters, biquadratic interac- 
tion, etc. It is worthwhile to remark that Minami [721471] has recently succeeded in obtaining a 
quite general exact solution for a rather large class of ID Ising models, which involve the most 
of aforelisted exactly solved models including the ones with the mixed spins and alternating 
bonds. 

2D Ising model had resisted almost two decades of intensive efforts until the complete closed- 
form exact solution has been found by Onsager for the spin-1/2 Ising model on a square lattice 
without the external magnetic field using the transfer- matrix method and Lie algebra [75J. 
Onsager's exact solution is currently regarded as one of the most significant achievements in the 
equilibrium statistical mechanics, because it brought an important revision in the understanding 
of phase transitions and critical phenomena. Actually, Onsager's exact results had served in 
evidence of a striking phase transition, which comes from extremely short-ranged interactions 
and is accompanied with a strange singular behaviour of several thermodynamic quantities in 
a close vicinity of the critical point. From this point of view, Onsager's exact solution also 
brought a considerable insight into deficiencies of some approximative methods, which mostly 
fail in predicting correct behaviour near a critical region. In this regard, one of the most 
essential questions to deal with in the statistical mechanics of exactly solvable models is always 
to find a precise nature of discontinuities and singularities accompanying each phase transition. 

The only disadvantage of Onsager's method lies in a considerable mathematical formida- 
bility of his solution. Bearing this in mind, many theoretical physicists have started to search 
for an alternative way, which would admit a more straightforward exact treatment of 2D Ising 
models. It is noteworthy that the original Onsager's solution has been slightly simplified by 
Kaufmann and Onsager himself using the theory of spinors [76|l77j. However, the more sub- 
stantial simplification has been later on achieved using various rigorous methods such as several 
combinatorial approaches developed by Kac, Ward and Potts [7j2[79], Hurst and Green [80J, 
Vdovichenko [8TT483] and others, the formalism of second quantization based on Jordan- Wigner 
fermionization invented by Schultz, Lieb and Mattis [84], the recurrence relations derived from 
the star-triangle transformation by Baxter and Enting [83], or more recent theories based on 
Grassmann variables [86-89] or Clifford-Dirac algebra [90] . Among all rigorous techniques, the 
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Pfaffian method seems to be the simplest method that enables to solve exactly the spin-1/2 
Ising model on 2D lattices. Within this rigorous approach, the problem of solving 2D Ising 
model is firstly reformulated as the problem of dimer statistics on a relevant decorated lat- 
tice [9Tj and subsequently, the relevant dimer statistics is precisely solved by employing the 
Pfaffian technique following the ideas of Kasteleyn [52], Temperley and Fisher [HHHM]. Note 
furthermore this exact method can be used for treating the spin-1/2 Ising model on any planar 
lattice without crossing bonds and among other matters, this method has thus proved a strong 
universality in a critical behaviour of 2D Ising lattices [95T498] . 

2.3 Heisenberg model 

It has been already mentioned in the preceding part that the quantum Heisenberg model 
[99J and its various extensions are much sought after, since they are more appropriate for 
modeling the magnetic behaviour of the most of real insulating magnetic materials than the 
aforedescribed Ising models. The main reason for a greater success of the Heisenberg-type 
models lies in the fact that these quantum spin models correctly take into account quantum 
fluctuations, which might play a crucial role in determining magnetic properties especially at 
low enough temperatures. In accordance with this statement, the magnetic behaviour of a 
rather extensive number of molecular-based magnetic compounds obeys theoretical predictions 
of the Heisenberg-type models as convincingly evidenced in the review article of de Jongh 
and Miedema [21], several earlier books devoted exclusively to the molecular-based magnetic 
materials |100f fl05j. as well as, the more recent series of books edited by Miller and Drillon 
[106f 4110j . It is safety to say that it is even not possible to enumerate manifold application of 
the Heisenberg-type models in explaining magnetic properties of real molecular-based magnetic 
compounds and hence, we will not dwell further on this aspect with respect to a rather extensive 
list of excellent literature [2"T }llUUIfllU] to be published on this topic. 

The quantum Heisenberg model can be defined through the Hamiltonian 

N 

n = -j^2s l -s J -HY / st, (2.2) 

where the first summation usually runs over all pairs of nearest-neighbour spins, the second 
summation is carried out over all lattice sites and the symbol Sj marks the spin operator for 
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the jth lattice site with the spatial components Sj = (Sj, Sj, Sj) to be given by two-by-two 
Pauli spin matrices^ 

S! = l\° 1 I , = i \ " | , s; = i \ ~ - | , (2.3) 



2 



1 





3 3 

where % = a/— T is the imaginary unit. As before, the former term in the Hamiltonian (12. 2p takes 
into consideration a configurational energy associated with the exchange interaction between 
the nearest-neighbour spins and the latter term accounts for the Zeeman's energy of discrete 
magnetic moments. In an attempt to treat exactly the quantum Heisenberg model one encoun- 
ters the most crucial mathematical difficulties in a non-commutability of the spin operators, 
which obey the following set of commutation relations 

[SJ, S k ] = S"S k — S^SJ = iSjSjkEap-y, (2.4) 

where a,/3,7 G {x, y, z}, the symbols 5jk and label Kronecker and Levi-Civita symbols, 
respectively. The most wide-spread extension of the quantum Heisenberg model (I2.2p certainly 
represents the XXZ Heisenberg model given by the Hamiltonian 

N 

U = -Jj2[ A (S?Sj + SfS]) + S-tS*] -Hj2Sf, (2.5) 

which merely replaces the isotropic pairwise spin-spin interaction through the anisotropic one. 
The anisotropy parameter A then allows to obtain the semi-classical Ising model as a special 
limiting case of the XXZ Heisenberg model in the limit A — > 0. Another special limiting case 
of the XXZ Heisenberg model, which is of particular research interest, represents the quantum 
XY model obtained from the Hamiltonian ( 12.51) by considering the other extreme limiting case 
A — >■ oo. 



2.4 Survey of exactly solved Heisenberg models 

Owing to mathematical complexities, which are closely connected with a non-commutability of 
the spin operators, the list of exactly solved quantum Heisenberg models is much less numerous 
compared to the above mentioned list of exactly solved Ising models. It is worthy of notice 
3 Here and in what follows the Pauli matrices are written in units % = 1. 
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that the exactly solved quantum Heisenberg models essentially comprise quantum spin models, 
which have been exactly solved only at zero temperature (ground state). Moreover, it should be 
also stressed that the quantum Heisenberg models exhibit a more diverse behaviour compared 
to the semi-classical Ising models just on behalf of the non-commutability of the spin operators, 
which in turn causes a presence of quantum fluctuations fundamentally influencing the ground 
state of the antiferromagnetic Heisenberg models. From this point of view, the ferromagnetic 
and antiferromagnetic quantum Heisenberg models differ basically one from another in their 
behaviour, the former exhibits a classical ferromagnetic ground state, whereas the latter one 
often exhibits a variety of unusual and exotic ground state(s) as the classical antiferromagnetic 
(Neel) order is not a true eigenstate of the Hamiltonians (12. 2p and (12. 5p . 

Bearing all this in mind, the antiferromagnetic quantum Heisenberg models remain at the 
forefront of theoretical research interest over the past few decades. The eigenstates and the 
ground-state energy of the antiferromagnetic spin-1/2 Heisenberg linear chain has been exactly 
found by Bethe [111] and Hulthen |112] with the help of the so-called Bethe-ansatz method. 
This method has been subsequently used also for a calculation of the magnetization [113] , 
the elementary excitation spectrum |114j . the susceptibility and other thermodynamic quan- 
tities [113] . It is noteworthy that the Bethe-ansatz method has been afterwards adapted to 
find similar exact solutions for the anisotropic spin-1/2 XXZ Heisenberg model on a linear 
chain as well |115H119] . Another interesting examples of the exactly solved quantum Heisen- 
berg models represent the antiferromagnetic spin-1/2 Heisenberg chains with the competing 
nearest-neighbour and the next-nearest-neighbour interactions, which are know as the so-called 
Majumdar- Ghosh model [120-123] and the delta-chain model |124H127] both having a rather 
peculiar dimerized ground state. Note furthermore that the rigorous approach, which was orig- 
inally developed by Majumdar and Ghosh [120J for the spin-1/2 Heisenberg chain with the 
competing nearest- and next-nearest-neighbour interactions, has been later on utilized when 
searching for an exact ground state of the antiferromagnetic spin-1/2 quantum Heisenberg 
model on several geometrically frustrated 2D and 3D lattices |128H133] . Last but not least, 
there also exist a rather large class of the exactly solved antiferromagnetic quantum Heisenberg 
models to be extended by various biquadratic and/or multispin interactions, which exhibit the 
intriguing ground state described within the valence-bond-solid picture [134H147] . 
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Finally, let us conclude our survey of the exactly solved quantum Heisenberg models by 
mentioning a few special limiting cases for which a quite general closed-form exact solution 
(not restricted only to the ground state) has been derived. Among these valuable exceptions 
one could mention the spin-1/2 quantum XY model on the linear chain |148H150] and the spin- 
1/2 Heisenberg-Ising bond alternating chain |148[I151] . which have been exactly treated using 
the Jordan- Wigner fermionization approach invented by Lieb, Schultz and Mattis [148J. Apart 
from these two fully exactly solvable quantum spin chains, the quite general exact analytical 
solution has been also found for the spin-1/2 quantum Heisenberg ladder with the intra-rung 
pair and inter-rung quartic interactions with the help of transfer- matrix method [152(1153] . 
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3 Dual transformation 

Some exact results for the spin-1/2 Ising model defined on 2D lattices can be obtained in a 
relatively straightforward way by making use of the dual transformation [51]. In this section, 
let us consider the spin-1/2 Ising model on 2D lattices in an absence of the term incorporating 
the effect of external magnetic field 

H = -JY J a i a j . (3.1) 

(hi) 

It is worthy to recall that cr, = ±1/2 represents two- valued Ising spin variable located at the 
ith lattice point and the summation is restricted only to all pairs of nearest-neighbour spins. 
Assuming that N is a total number of lattice sites and z is being its coordination number 
(i.e. the number of nearest neighbours), then, there is in total = Nz/2 pairs of nearest- 
neighbour spins when boundary effects are neglected in the thermodynamic limit iV — > oo. Each 
line (bond), which connects two adjacent spins on 2D lattice, can be regarded as a schematic 
representation of the pairwise interaction J between the nearest-neighbour spins. As usual, the 
central issue of our approach is to calculate the configurational partition function 

Z=$>xp(-/3ft), (3.2) 
R} 

where /3 = 1/(/cbT), &b is Boltzmann's constant, T is the absolute temperature and the suffix 
{cTj} denotes a summation over all possible configurations of the Ising spins on a given 2D 
lattice. First, let us rewrite Hamiltonian (13.1 j) to the form 




Since the canonical ensemble average of Hamiltonian readily represents the internal energy, 
i.e. U = (H), it is easy to find the following physical interpretation of the Hamiltonian ( 13. 3ft . 
Each couple of unlike oriented adjacent spins contributes to the sum on the right-hand-side of 
Eq. (13.31) by the energy gain J / 2, while each couple of aligned adjacent spins does not contribute 
to this sum at all. With respect to this, the internal energy of a completely ordered spin system 
(all spins either 'up' or 'down') acquires its minimum value U = —N-bJ/4 under the assumption 
J > 0. By substituting the Hamiltonian (13.30 into Eq. (13.21) . one gets the following expression 
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for the partition function 



Z = exp ^— 




) 



(3.4) 



where n stands for the total number of unaligned spin pairs within each spin configuration. 
Obviously, the sum on the right- hand- side of Eq. H3.4[) may in principle contain just different 
powers of the expression exp (—(3 J/2). In the limit of zero temperature (T — > 0), the expression 
exp(— j3J/2) tends to zero and hence, the power expansion into a series ex P( — nj3J/2) gives 



very valuable estimate of the partition function in the limit of sufficiently low temperatures, 
the low-temperature series expansion [51~ |ll54l 4l56j. 

It is possible to find a simple geometric interpretation of each term emerging in the afore- 
mentioned power series by introducing a dual lattice. For this purpose, let us briefly mention a 
basic terminology of the graph theory [157( 1158] . where each site of a lattice is called as a vertex, 
while each bond (line) connecting the nearest-neighbour sites (vertices) is called as an edge. 
Further, an interior of each elementary polygon delimited by edges is a face and an ensemble of 
vertices and edges is called a full lattice graph. Vertices of a dual lattice are simply obtained by 
situating them in the middle of each face of the original lattice. The vertices situated at adjacent 
faces, which share a common edge on the original lattice, then represent the nearest-neighbour 
vertices of the dual lattice. Edges of the dual lattice are obtained by connecting each couple of 
adjacent vertices of the dual lattice. For illustrative purposes, Fig. [1] shows two original lattices 
- square and honeycomb - and their corresponding dual lattices. The bonds of original lattices 
are displayed in this figure as solid lines, while the bonds of their dual lattices are depicted as 
broken lines. As one can see, the square lattice is a self-dual, i.e., the dual lattice to a square 
lattice is again a square lattice. Contrary to this, the triangular lattice is a dual lattice to the 
honeycomb lattice and vice versa. Remembering that N and is the total number of vertices 
and edges of the original lattice, respectively, it follows that is at the same time the total 
number of edges of the dual lattice as well (each bond of the dual lattice intersects one and 
just one bond of the original lattice). If Nn denotes the total number of vertices of the dual 
lattice, then, Nd also equals to the total number of faces of the original lattice (each face of the 
original lattice involves one and just one vertex of the dual lattice). Euler's relation for planar 



n 
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a) b) 

Figure 1: Two pairs of dual lattices: a) self-dual square lattices; b) honeycomb and triangular dual 
lattices. Solid lines and solid circles label edges and vertices of original lattices, while broken lines and 
empty circles stand for edges and vertices of their dual lattices, respectively. 

graphgj consequently relates the total number of vertices of the original and dual lattices with 
the total number of edges 



N + N D = N E 



(3.5) 



Suppose now a random arrangement of 'up' and 'down' spins on the original lattice as it 
is displayed in Fig. [2] for some particular example of a spin configuration on a square lattice. 
If the pairs of adjacent spins are not aligned alike, then draw solid lines between them, other- 
wise draw broken lines between them. The system of vertices, solid and broken lines creates 
a configurational graph on the dual lattice, which is a subgraph of the full dual lattice graph. 
The most fundamental property of the configurational graph is that the mutual interchange of 
each couple of adjacent spins either does not affect the configurational graph, or it causes an 
even number of changes. The total number of solid and broken lines incident at each site of the 
dual lattice must be therefore either even or zero. This implies that solid (broken) lines of each 
configurational graph form a system of closed polygons and hence, it follows that the configura- 
tional graph represents a certain kind of polygon line graph. Another important property of the 
polygon line graph is that the reversal of all spins does not change the configurational graph. 
It means that two different spin configurations (one is obtained from the other by reversing all 



The planar graph is a graph, which can be embedded in the plane so that no edges intersect. 
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Figure 2: A particular spin configuration on a square lattice, whose vertices and edges are not drawn 
for clarity. The plus (minus) sign at ith lattice site corresponds to the spin state cr, = +1/2 (—1/2). 
The system of solid and broken lines unambiguously determines the corresponding polygon line graph 
on the dual square lattice (for details see the text). 

spins) correspond just to one configurational polygon line graph due to invariance cr, — > — (Tj 
(i = 1, 2, ... , N). The partition function of 2D Ising model can be therefore expressed as 

(3J\ 



Z = 2exp ^iVB^exp 



-n- 



(3.6) 



where the summation is now performed over all possible polygon subgraphs on the dual lattice, 
n denotes the total number of solid lines within each polygon subgraph and the factor 2 comes 
from the two-to-one mapping between spin and polygon configurations. 

Now, let us take a closer look at another interesting property of the partition function (13 .2p . 
By adopting the exact van der Waerden identity [159J 



exp(/5 Jai<jj 



cosh f~j~ j + 4cTjOj sinh f 



cosh 



1 + 4:<Ji<jj tanh I — ^- 



(3.7) 



and substituting it into Eq. (13. 2 h one obtains 



cosh 



/3J 



En 



1 + 4<7j<7j tanh 



(3.8) 



The product on the right-hand-side of Eq. (13. 8 \ involves in total N-q terms, which give after 
formal multiplication a sum of in total 2^ terms. However, many terms eventually vanish after 
performing a summation over all available spin configurations. For instance, it can be readily 
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proved that all linear terms of the type OjOj tanh(/3J/4) will disappear after summing over spin 
states of either the spin Oj or o~j. The condition, which ensures that the relevant term makes a 
non-zero contribution to the partition function, can be simply guessed from the validity of the 
trivial identity of = 1/4. Accordingly, all non-zero terms must necessarily contain only spin 
variables, which enter into these expressions either even number of times or do not enter into 
these terms at all. The simplest non-vanishing term for the Ising square lattice is evidently the 
term 



l 4 (oWj) (0jcr fc ) (a k ai) (aiai) tanh 4 f-r- j = 4 4 of o| of of tanh 4 i-^- j = tanh 4 



which is constituted by the product of four nearest-neighbour interactions (bonds) whose cor- 
responding edges form an elementary square (the simplest closed polygon) on this lattice. The 
summation over spin configurations of four spins included in this term consequently gives the 
factor 2 4 tanh 4 ((3 J/4), while the summation over spin states of other spins yields the additional 
factor 2 N ~ 4 so that the Boltzmann factor 2 N tanh 4 (/3J/4) is finally obtained as the contribution 
from a single square. It is therefore not difficult to construct the following geometric interpre- 
tation of the non-vanishing terms: edges corresponding to interactions to be present in these 
terms must create closed polygons so that either no lines or even number of lines meet at each 
vertex of 2D lattice. From this point of view, polygon line graphs very similar to those described 
by the low-temperature series expansion give non-zero contributions to the partition function. 
With regard to this, the expression (13. 8p for the partition function can be replaced with 



Z = 2 N 



cosh 



Nb r //?j x 

V] tanh 
p-g- 



(3.9) 



where n denotes the total number of full lines constituting the particular polygon graph. It is 
quite apparent that the summation on the right-hand-side of Eq. (13. 9 j) contains just different 
powers of the expression tanh(/3 J/4). In the limit of high temperatures (T — > oo), the expression 
tanh(/3J/4) tends to zero and hence, the power expansion into a series ^[tanh(/3 J/4)] 71 gives 

n 

very valuable estimate of the partition function in the limit of high enough temperatures, the 
high-temperature series expansion [511115414156] . 

There is an interesting correspondence between summations to emerge in Eqs. (13.61) and 
(13. 9p . since both of them are performed over certain sets of polygon line graphs. The most 
essential difference between them lies in the fact that the summation in Eq. (13 .4p is carried out 
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over polygon graphs on the dual lattice, while the summation in Eq. ( 13.91) is performed over 
polygon graphs on the original lattice. It should be stressed, however, that both expressions 
(13. 4p and (13.91) for the partition function are exact when the relevant series is performed up 
to an infinite order and therefore, they must basically give the same result for the partition 
function. In the thermodynamic limit, the summations (13 .4p and (13. 9p yield the same partition 
function provided that 



cxp 



for 



Z(N D ,f3 D J) 



Z(N,f3J) 



exp 



cosh (¥) 



(3.10) 
(3.H) 



where we have introduced the reciprocal temperature fo = l/(k B T D ) of the dual lattice and 
the factor 2 was omitted from the denominator on the left-hand-side of Eq. (13. lip as it can 
be neglected in the thermodynamic limit. The connection between two mutually dual lattices 
( 13.101) and ( I3.1ip can also be inverted because of a symmetry in the duality 



cxp 
Z{N,PJ) 



¥ — 



Z(N D ,foJ) 



exp (^iV B ) 2^ [cosh (^) 



(3.12) 
(3.13) 



since each from a couple of the mutually dual lattices is dual one to another. With the help 
of Eqs. (I3.10p and (13. lip [or equivalently Eqs. (I3.12p and (I3.13p ]. it is also possible to write 
this so-called dual transformation even in a symmetric form as it could be expected from the 
symmetrical nature of the duality. The relation ( I3.10p gives after straightforward algebraic 
manipulation 



sinn — smh 

\ 2 \ 2 . 



(3.14) 



while the expression (13. lip can be modified by regarding the equality between partition func- 
tions, Eqs. ( 13.1 Op . (I3.14p and Euler's relation (I3.5P implying that 



cosh 



cxp 



Nv 



2sinh(^) 



_> 



(3.15) 



V 1 [2sinh(^y 

By combining Eq. (13.151) with the relation (13.111) . one actually gains another symmetric rela- 
tionship between the partition functions Z(N,(3J) and Z(Nd, foJ), which are expressed in 
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terms of the high-temperature series expansion on the original lattice and the low-temperature 
series expansion on its dual lattice 

Z{N,/3J) Z(N D ,(3 D J) 



N 



(3.16) 



2sinh(f)] 2 [2sinh(^)" 

An existence of the mapping equivalence between the low- and high-temperature series expan- 
sions reflects the fundamental property of the partition function, namely, its symmetry with 
respect to the low and high temperatures. This symmetry means, among other matters, that 
the partition function at some lower temperature can always be mapped on the equivalent 
partition function at some certain higher temperature. This mapping is called as the dual 
transformation and the dual lattices are actually connected one to another by means of the 
dual transformation. In this respect, the dual lattices are topological representations of the 
dual transformation and consequently, one says that the dual transformation has a character 
of the topological transformation. 

The mathematical formulation of the dual transformation connecting effective temperatures 
of the original and its dual lattice is represented (independently of the lattice topology) either 
by the couple of equivalent equations ( 13 . 1 j) and (I3.12p . or, respectively, by a single symmetrized 
relation ( I3.14p . The latter relationship is especially useful for a better understanding of the 
symmetry of the partition function with respect to the low and high temperatures. It is suffi- 
cient to realize that the argument of the function sinh(/3J/2) must unavoidably decrease when 
the argument of the other function sinh(/3D<//2) increases in order to preserve their constant 
product required by the dual transformation (I3.14p . This means that the partition function at 
some lower temperature, which is obtained for instance from the low-temperature expansion 
on the dual lattice, is equivalent to the partition function at a certain higher temperature ob- 
tained from the high-temperature expansion on the original lattice. It is noteworthy that the 
dual transformation markedly simplifies the exact enumeration of thermodynamic quantities on 
different lattices, since it permits to obtain the exact solution of some quantity on arbitrary 2D 
lattice merely from the corresponding exact result for this quantity on its dual 2D lattice. For 
instance, the critical point is always accompanied with some singularity or discontinuity in the 
thermodynamic functions and this non-analyticity is always somehow reflected in the partition 
function as well. The mapping relation (I3.16P clearly shows that the partition function of the 

18 



3 DUAL TRANSFORMATION 



original lattice exhibits a non-analyticity if and only if the partition function of the dual lattice 
also has a similar non-analyticity at some corresponding temperature satisfying the duality 
relation (I3.14p . Besides, the expression (I3.16P allows to calculate the partition function of the 
spin-1/2 Ising model on the one from two mutually dual lattices merely from the corresponding 
exact result for the partition function of its dual lattice model. 

It is worthwhile to remark that the square lattice has an extraordinary position among 2D 
lattices because of its self-duality. The self-dual property together with the symmetry of the 
partition function with respect to the low and high temperatures is just enough for determin- 
ing the critical temperature and other thermodynamic quantities precisely at a critical point. 
Under the assumption of a single critical point, the same lattice topology ensures that critical 
parameters must be equal one to another on both mutually dual square lattices. According 
to the dual transformation (13.141) . the critical point of the spin-1/2 Ising model on the square 
lattice must obey the condition 

sinh 2 (^fj = 1, (3.17) 

which is consistent with this value of the critical temperature T c [f3 c = l/(k B T c )} 
kuT c 1 

-ttt = 7— ■ 3.18 

\J\ 2 ln(l + >/2) 

The location of a critical point of the spin-1/2 Ising model on the square lattice, which has 
been achieved in 1941 by Kramers and Wannier [51 J with the help of the dual transformation, 
can be regarded as the first exact analytical result serving in evidence of the spontaneous long- 
range ordering. However, the main disadvantage of the dual transformation lies in the fact that 
this exact method cannot serve for calculating thermodynamic quantities out of the critical 
point and moreover, the temperature symmetry of the partition function does not suffice for 
determining critical parameters of other 2D lattices like honeycomb and triangular lattices, 
which are not self-dual. 
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4 Algebraic transformations 

It has been demonstrated in the preceding section that it is impossible to find from the dual 
transformation alone an exact critical point of the spin-1/2 Ising model on 2D lattices, which 
are not self-dual. This has stimulated considerable interest in the search for algebraic mapping 
transformations, which would tackle this outstanding problem when combining them with the 
dual transformation. In what follows, our attention will be therefore focused on the most 
important features of algebraic transformations. 

4.1 Star-triangle transformation 

The dual transformation (l3.14j) maps the honeycomb lattice into the triangular lattice or vice 
versa and thus, it does not establish the symmetry between the low- and high-temperature 
partition function on the same lattice. It is worthy of notice that such a precise relation can 
be alternatively derived by making use of some algebraic mapping transformation. The star- 
triangle transformation, which was originally invented by Onsager in his famous work |75j, 
establishes this useful mapping relationship for two interesting couples of dual lattices such as 
honeycomb-triangular and kagome-diced lattices. 

Consider the spin-1/2 Ising model on the honeycomb lattice of 2N sites. For further con- 
venience, it is advisable to divide the honeycomb lattice into two equivalent interpenetrating 
triangular sublattices, whose sites are diagrammatically represented in Fig. [3] by full and empty 
circles, respectively. The division is made in a such way that all nearest neighbours of a site 
from the first sublattice belong to the second sub lattice and vice versa. Owing to this fact, the 
summation over configurations of spins, which belong to the same sublattice, can be performed 
independently one from each other because of absence of any direct interaction between the 
spins from the same sublattice. For easy reference, let us formally denote the spins from the 
first sublattice as the vertex spins o~i and the spins from the second sublattice as the decorat- 
ing spins Hi. The total Hamiltonian can be for further convenience written as a sum of site 
Hamiltonians 

N 

1=1 
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4.1 Star-triangle transformation 




STT 



• a, 



Figure 3: The spin- 1/2 Ising model on the honeycomb lattice and its relation to the equivalent spin- 
1/2 Ising model on the triangular lattice. The equivalence between both models can be established 
by applying the star-triangle transformation to a half of spins of honeycomb lattice. 

where each particular site Hamiltonian TLi involves all the interaction terms associated with 
the ith decorating spin 



Hi = —Jh^iivn + o"i2 + cr 



/3J 



(4.2) 



By the use of Eqs. (14.11) and (14. 2p . the partition function of the spin- 1/2 Ising model on the 
honeycomb lattice (13.21) can be partially factorized to the form 



N 



Z ^ = Y II Y exp[/3J h /!i(cTa + a i2 + cr i3 )], 

{<Ti} 1=1 tn=±l/2 



(4.3) 



where the former summation is carried out over all available configurations of the vertex spins, 
the product runs over all decorating spins and the latter summation accounts for the spin 
states of one particular decorating spin. Consequently, it is adequate to consider the individual 
decorating spin /i; and to sum up over degrees of freedom of this spin. It should be stressed 
that each decorating spin interacts merely with its three nearest-neighbour vertex spins and 
hence, this summation gives the Boltzmann's factor 

~/3J h . 



exp\pj h fii(an + a i2 + a i3 )] = 2 cosh 

W=±l/2 



0"ii + &i2 + cr. 



Aexp[/3Jt(a il ai2 + cr i2 a i3 + a^an)], 



(4.4) 



which can be substituted by a simpler equivalent expression provided by the star-triangle trans- 
formation. The physical meaning of the mapping transformation (14.41) lies in removing all the 
interaction terms associated with the central decorating spin m of the star and replacing them 
with some effective interactions between the three outer vertex spins an, 012 and forming 
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the equilateral triangle. It is noteworthy that the star-triangle transformation is actually a set 
of eight equations, which can be obtained from the algebraic transformation (14. 4p by consid- 
ering all possible spin configurations available to the three outer vertex spins. However, the 
consideration of eight available spin configurations leads just to two independent equations, 
which unambiguously determine the mapping parameters A and f3J t 



A 



0Jt = ^ 



cosh 



cosh f 



4 J 

3P A \ 



cosh 



'PJh 



cosh 



(4.5) 



(4.6) 



The backward substitution of the transformation (|4.4j) into the partition function (14. 3p . which 
is equivalent with performing the star-triangle transformation for all decorating spins /ij, yields 
an exact mapping relationship between the partition functions of the spin- 1/2 Ising model on 
the honeycomb and triangular lattices 

Z h (2N, 0J h ) = A N Z t (N, 0J t ), (4.7) 

whose corresponding temperatures are coupled together through the mapping relation (14. 6p . 
As a result, the mapping relation ( 14. 6p resulting from the star-triangle transformation connects 
the partition functions of the honeycomb and triangular lattices at two different temperatures 
in a very similar way as it does the relation (13. 14|) provided by the dual transformation. The 
most crucial difference consists in a profound essence of both mapping relations; the dual 
transformation is evidently of the topological origin, whereas the star-triangle mapping is the 
algebraic transformation in its character. 

At this stage, let us combine the dual and star-triangle transformations in order to bring 
insight into a criticality of the spin- 1/2 Ising model on the honeycomb and triangular lattices. 
By employing a set of trivial identities for hyperbolic functions, the star-triangle transformation 
(14.61) can also be rewritten as follows 



cosh f (BJC 
exp (/3J t ) = = 2 cosh ' h 



- 1. 



(4.8) 



cosh(^) V 2 , 

Furthermore, it is appropriate to combine Eq. (14. 8 j) with one of possible representations of the 
dual transformation (I3.10p 



exp 



P'Jt 



tanh 



(4.9) 
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4.1 Star-triangle transformation 



with the aim to eliminate temperature of the one from two lattices with the equivalent partition 
functions. For instance, the procedure that eliminates from Eqs. (14. 8p and (14. 9p the effective 
temperature of the triangular lattice allows one to obtain the symmetrized relationship, which 
connects the partition function of the spin-1/2 Ising model on the honeycomb lattice at two 
different temperatures 



cosh 



cosh 



'Ft 



ii 



1. (4.10) 



Note that the relationship (I4.10p establishes analogous temperature symmetry for the partition 
function of the spin-1/2 Ising model on the honeycomb lattice as the dual transformation does 
for the spin-1/2 Ising model on the self-dual square lattice through the relation (I3.14p . If there 
exists just an unique critical point, both temperatures connected via the mapping relation 
(14.101) must necessarily meet at a critical point due to the same reason as it has already been 
explained by analysis of the square lattice. So, the critical point of the spin-1/2 Ising model on 
the honeycomb lattice should obey the condition 



(4.11) 



(4.12) 



COsh ( j - 1 

which is consistent with this value of the critical temperature 

kyT c _ 1 
|J h | ~ 21n(2 + y/3)' 

The same procedure can be repeated once more in order to obtain the critical parameters of 
the spin-1/2 Ising model on the triangular lattice. An elimination of the effective temperature 
of the honeycomb lattice from Eqs. ( 14. 8ft and ( 14.91) yields the following symmetric relationship 

[exp(/3J t ) - l][exp(/3'J t ) - 1] = 4, (4.13) 

which relates the partition function of the spin-1/2 Ising model on the triangular lattice at two 
different temperatures. The relation (14. 13ft consecutively determines the critical condition for 
the spin-1/2 Ising model on the triangular lattice 

[exp(/3 c J t ) - l] 2 = 4, (4.14) 



which locates its exact critical temperature 

fcB^c _ _J_ 

J t ~ In 3' 
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Figure 4: The spin-1/2 Ising model on the honeycomb lattice, decorated honeycomb lattice and 
kagome lattice. The equivalence between all three lattice models can be established by employing the 
decoration-iteration (DIT) and star-triangle (STT) transformations, respectively. 

It is noteworthy that the critical temperature of the spin-1/2 Ising model on the triangular 
lattice (I4.15P can be more easily found by substituting the critical condition of the spin-1/2 
Ising model on the honeycomb lattice (I4.1ip to the star-triangle transformation ( I4.8p . Of course, 
the same critical temperature will be recovered also in this way. 

4.2 Decoration-iteration transformation 

Another important mapping transformation, which is purely of algebraic character, represents 
the decoration-iteration transformation invented by Syozi in 1951 by solving the spin-1/2 Ising 
model on the kagome lattice |160j . In principle, the approach based on the decoration-iteration 
transformation enables to obtain an exact solution of the spin-1/2 Ising model on an arbitrary 
bond decorated lattice from the corresponding exact solution of the simple (undecorated) lattice. 
The term bond decorated lattice marks such a lattice, which can be obtained from a simple 
original lattice (like square, honeycomb, triangular or any other) by placing an additional spin 
or a finite cluster of spins on the bonds of this simpler lattice. The spins placed at vertices of 
the original lattice will be referred to as vertex spins, while the additional spins arising from the 
decoration procedure will be further called as decorating spins. For illustration, Fig. H] shows 
the planar topology of the honeycomb lattice, the simply decorated honeycomb lattice and the 
kagome lattice together with the mapping relations, which can be established between them by 
employing algebraic mapping transformations. 

Let us consider initially the spin-1/2 Ising model on the decorated honeycomb lattice, which 
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4.2 Decoration-iteration transformation 



is shown in the middle of Fig. H] and is given by the Hamiltonian 

3N 

n = -J d J2^y (4-16) 

(v) 

Above, Oi = ±1/2 and fij = ±1/2 label the vertex and decorating Ising spins, respectively, 
the summation runs over all nearest-neighbour spin pairs on the decorated honeycomb lattice 
and the total number of the vertex spins is set to N. It is convenient to rewrite the total 
Hamiltonian (I4.16P as a sum of bond Hamiltonians 

3JV/2 

n = E % (4-i7) 

3=1 

where each bond Hamiltonian Hj involves all the interaction terms of the decorating spin fij 
from the jth bond of the decorated honeycomb lattice 

Hj = -J d /ij-(0-ji + cr j2 ). (4.18) 

By the use of Eqs. (14.171) and (14.181) . the partition function (13.21) of the spin-1/2 Ising model 
on the decorated honeycomb lattice can be partially factorized to the form 

37V/2 

^d = E II E exp\PJdl*j{°fl. + (4.19) 

Wi} j=l /Uj=±l/2 

In the above expression, the former summation is carried out over all available configurations of 
the vertex spins, the product runs over all decorating spins and the latter summation accounts 
for spin states of the decorating spin rij. According to Eq. ( I4.19p . the summation over degrees 
of freedom of the decorating spins can be performed independently one from each other and 
hence, this summation gives the effective Boltzmann's factor 



J2 exp^Jd/x^cTji + cr j2 )] = 2 cosh 



{(Tjl + CT j2 ) 



BexptfJ&jKTv), (4.20) 



which can be successively replaced with a simpler equivalent expression provided by the decoration- 
iteration transformation. The physical meaning of the mapping transformation (14.201) lies in 
removing all the interaction terms associated with the decorating spin rij and substituting them 
by the effective interaction between the two vertex spins o~j\ and aj 2 , which are being its near- 
est neighbours. It should be emphasized that the decoration-iteration transformation (14.201) 
has to satisfy the self-consistency condition, i.e., it must hold for any combination of the spin 
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states of the two vertex Ising spins <jj\ and 0^2 • In this respect, the mapping relation (I4.20p 
is in fact a set of four equations, which can be explicitly obtained by considering all possible 
spin configurations available to the two vertex spins Oj\ and Oj%. It can be readily proved that 
a substitution of four available spin configurations yields from the formula (I4.20p merely two 
independent equations, which determine so far not specified mapping parameters B and fiJ^ 



B 



cosh 



/3 J h = 2 In 



cosh 



2 J 



(4.21) 
(4.22) 



2 

By applying the decoration-iteration transformation to all decorating spins, i.e. substituting 
the mapping transformation (I4.20p into the partition function (14.191) . one acquires an exact 
mapping correspondence between the partition functions of the spin- 1/2 Ising model on the 
decorated honeycomb lattice and simple honeycomb lattice 

Z d (5N/2, (3J dec ) = B m / 2 Z h (N, 0J h ), (4.23) 

whose effective temperatures are connected by means of the relation (14.221) . It is quite obvi- 
ous from Eq. (14.231) that the decorated honeycomb lattice becomes critical if and only if its 
corresponding honeycomb lattice becomes critical as well. With regard to this, it is sufficient 
to substitute the exact critical temperature of the honeycomb lattice f 14 . 1 2 1) into Eq. (I4.22[) in 
order to locate the critical point of the decorated honeycomb lattice 

knT c 1 

ITT = 7- — / ( 424 ) 

l J dl 2 ln(2 + V3 + ^6 + 4^3) 
It is worthwhile to remark that the decoration-iteration transformation ( 14 . 2 [) is not restricted 
neither by the geometry of a lattice nor by its spatial dimensionality and thus, it can be 
utilized for obtaining rigorous results for arbitrary simply decorated lattice from the known 
exact solution of the corresponding undecorated lattice. 

Now, it is possible to make another useful observation. The total Hamiltonian (14. 16j) of the 
spin-1/2 Ising model on the decorated honeycomb lattice can also be formally written as a sum 
of site Hamiltonians 

N 

« = ( 4 - 25 ) 

1=1 
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4.2 Decoration-iteration transformation 



whereas each particular site Hamiltonian Hi involves all the interaction terms of the one indi- 
vidual vertex Ising spin <7j 



/3J 



(4.26) 



Substituting Eqs. (14.251) and (I4.26P into a statistical definition of the partition function (13.21) 
allows one to partially factorize the partition function of the spin- 1/2 Ising model on the 
decorated honeycomb lattice and to rewrite it into the form 



N 



Z a = II Yl exp\pJi<Ti(Pii + A*i2 + Ha)]- 

{fli} 1=1 <7i=±l/2 



(4.27) 



Above, the former summation runs over all available configurations of the decorating spins, 
the product runs over all vertex spins and the latter summation is carried out over the spin 
states of the ith vertex spin cr, from the decorated honeycomb lattice. The structure of the 
partition function (I4.27P immediately justifies applicability of the familiar star-triangle mapping 
transformation 

~0Jd, 



exp\pj A <Ti{pa + Hi 2 + /Ms)] = 2 cosh 

(Ti=±l/2 



■(Hn + u, i2 + A* 



a) 



CexTp[/3J k (u.a/ii2 + ^2^3 + ^3^1)], 



which satisfies the self-consistency condition provided that 



C 



fij k = In 



cosh 



'3/3 J/ 



cosh (2^1 j 
cosh 



cosh 



In 



2 cosh 



'PJd 



(4.28) 

(4.29) 
(4.30) 



The star-triangle transformation maps the spin- 1/2 Ising model on the decorated honeycomb 
lattice into the spin- 1/2 Ising model on the kagome lattice once it is performed for all vertex 
spins (7j. As a matter of fact, it is easy to derive the following exact mapping equivalence 
between the partition functions of both these models by a direct substitution of the star-triangle 
transformation (14.281) into the relation (I4.27P 



Z d (5iV/2,/3J d ) = C N Z k (3N/2,l3J k ). 



(4.31) 



At this stage, it is possible to combine the decoration-iteration transformation with the 
star-triangle transformation in order to express the partition function of the spin- 1/2 Ising 
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honeycomb 


kagome 


square 


triangular 


k B T c / J 


0.37966 


0.53583 


0.56729 


0.91024 



Table 1: Critical temperatures of the spin- 1/2 Ising model on several 2D lattices. 

model on the kagome lattice through the corresponding partition function of the spin-1/2 Ising 
model on the honeycomb lattice. The mapping relations (I4.23P and (I4.3ip provide this useful 
connection between the partition functions of the spin-1/2 Ising model on the honeycomb and 
kagome lattices 



Z k (3N/2,f3J k ) = 2 
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exp (Mi) 
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2 exp(£i 
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exp 
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Z h (N,pj h ), 



(4.32) 



whereas Eqs. (14.221) and (14.301) relate the effective temperatures of the honeycomb and kagome 
lattices with the equivalent partition functions 



/3J k = ln 



2 exp 



(4.33) 



Substituting the exact critical temperature of the spin-1/2 Ising model on the honeycomb lattice 
(14.121) to the mapping relation (I4.33j) yields the exact critical temperature of the spin-1/2 Ising 



model on the kagome lattice 
k B T c 1 



(4.34) 



4 ln(3 + 2\/3)' 

Before proceeding further, let us compare the obtained critical temperatures corresponding 
to the order-disorder phase transition of the spin-1/2 Ising model on several 2D lattices. For this 
purpose, the Table [T] enumerates critical temperatures of the spin-1/2 Ising model on all three 
regular planar lattices - honeycomb, square and triangular, which are the only particular plane 
tilings that entirely cover the whole plane with the same regular polygon - hexagon, square 
and triangle, respectively [161] . The critical point of the semi-regular kagome lattice, which 
consists of two kinds of regularly alternating polygons (hexagons and triangles) is especially 
interesting from the academic point of view, since this 2D lattice represents the only semi- 
regular tiling, which has all sites as well as all bonds equivalent quite similarly as a triad of the 
aforementioned regular lattices. It can be easily understood from the Table [T] that the greater 
the coordination number (the number of nearest neighbours) of the planar lattice is, the higher 
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is the critical temperature of its order- disorder transition. Accordingly, the cooperativity of 
spontaneous ordering seems to be very closely connected with such a topological feature as 
the coordination number of 2D lattice is. On the other hand, the coordination number by 
itself does not entirely determine the critical temperature as it can be clearly seen from a 
comparison of the critical temperatures of the square and kagome lattices having the same 
coordination number four, but slightly different critical temperatures. It turns out that the 
semi-regular Archimedean lattices |161j composed of more regular polygons always have lower 
critical temperature than their regular counterparts with the same coordination number. Thus, 
it might be concluded that the cooperativity is strongly related also to other topological features 
of planar lattices. It should be also noted here that such an information cannot be elucidated 
from rough approximative methods, because the most of them usually predict the same critical 
temperature for 2D lattices with the same coordination number. 

4.3 Generalized transformations I 

In two preceding parts of this section, we have shown usefulness of algebraic mapping trans- 
formations in providing several exact results for the spin- 1/2 Ising model on 2D lattices after 
performing relatively modest calculations. From this point of view, the quite natural ques- 
tion arises whether or not algebraic transformations can be further extended and generalized. 
It is worth noticing that an early development of the concept based on generalized algebraic 
transformations has been elaborated by Fisher in the comprehensive paper to be published 
more than a half century ago [162] . In this work, Fisher has questioned a possibility of how 
the decoration-iteration transformation, the star-triangle transformation and other algebraic 
transformations can be generalized and besides, this notable paper has also furnished a rigor- 
ous proof on a general validity of algebraic mapping transformations. It should be remembered 
that algebraic transformations are carried out at the level of the partition function and their 
physical meaning lies in replacing a conveniently chosen part of the partition function with 
a simpler equivalent expression (Boltzmann's factor) to be obtained after performing a trace 
over degrees of freedom of a single decorating spin or a finite number of decorating spins, re- 
spectively. It is of principal importance that summations over spin states of the decorating 
spins (or the finite cluster of decorating spins) can be performed independently one from each 
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Figure 5: A diagrammatic representation of the generalized decoration-iteration transformation (DIT) 
for the decorating system composed of a single decorating spin S k of arbitrary magnitude. The terms 
Ji and J 2 stand for the pair interactions between the decorating spin S k and its two nearest-neighbour 
vertex spins cr k i arid a k2 , while the terms Hs and H a represent the effect of external magnetic field 
acting on the decorating and vertex spins, respectively. 

other and before summing over spin states of the vertex Ising spins. As a result, the validity 
of generalized algebraic transformations can be verified locally by considering the Boltzmann's 
factor of the relevant spin cluster, which consists of the central decorating spin (or the finite 
number of decorating spins) coupled to a few outer vertex Ising spins. 

4.3.1 Generalized decoration-iteration transformation 

Consider first the generalization of the decoration-iteration transformation, which is applicable 
for the models in which each decorating system interacts merely with the two outer vertex Ising 
spins. Let a single decorating Ising spin Sk of arbitrary magnitude be the decorating system, 
which interacts with the two vertex spins a k \ and a k2 as it is schematically illustrated in Fig. |5j 
Note that the subsequent generalization to a more general decorating system composed of a 
finite number of decorating spins is quite straightforward and will be investigated hereafter. The 
generalized decoration-iteration transformation for the decorating system, which constitutes a 
single decorating Ising spin Sk of arbitrary magnitude, can be mathematically formulated as 
follows 

s 

exp[(3H a (a kl + a k2 )] ^ exp [f3S k {Jia kl + J 2 cr k2 + H s )} 

S k =—S 

= Aexp[f3Ra kl cr k2 + 0H 1( T kl + f3H 2 a k2 \. (4.35) 

It is worthwhile to remark that the expression on the left-hand-side of the generalized decoration- 
iteration transformation ( I4.35P is in fact the Boltzmann's factor of the three-spin cluster from 
the left-hand-side of Fig. \5\ which enters into the partition function of the Ising model tak- 
ing into account two different pair interactions J\ and J 2 between the central decorating spin 
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and two outer vertex spins, as well as, the magnetostatic Zeeman's energy of the decorating 
and vertex spins in a presence of the external magnetic field of the magnitude Hs and H a , 
respectively. This expression is substituted through the generalized mapping transformation 
( I4.35P by the multiplicative factor A and the Boltzmann's factor of the two-spin cluster from 
the right-hand-side of Fig. which takes into account the pair interaction R between the two 
outer vertex spins, as well as, the effect of generally non-uniform magnetic field Hi and H 2 
acting on the vertex spins <Tki and cr^, respectively. The physical meaning of the decoration- 
iteration transformation (I4.35P lies in removing all the interaction terms involving the central 
decorating spin Sk and substituting them through a simpler equivalent expression depending 
just on the two vertex Ising spins. It is noteworthy that the generalized decoration- iteration 
transformation holds quite generally if and only if the mapping transformation (I4.35P remains 
valid regardless of four different spin configurations of the two vertex spins er^i and a^ 2 . This 
self-consistency condition determines so far not specified mapping parameters 



a = {y^v^Y 



PR = In 



V 3 V 4 



1 /ViVa 
^ 2 = + I In (|| 



(4.36) 



which are expressed in terms of the functions Vj (j 
abbreviated and elegant form 



1 — 4) in order to write them in a more 
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(4.37) 



It should be stressed that the functions Vj (j 



4) are actually four different Boltzmann's 



weights, which can be obtained from the Boltzmann's factor on the left-hand-side of the general- 
ized decoration-iteration transformation (14.351) by considering four possible spin configurations 
available to the two vertex Ising spins a^i and a^. 

At this stage, let us consider a more general decorating system composed of a finite number 
of the decorating spins instead of a single decorating spin as it is schematically depicted in 
Fig. EJ It has already been pointed out by Fisher [162J that a single decorating spin can be 
in principle replaced with an arbitrary (even quantum) statistical-mechanical system without 
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Figure 6: A diagrammatic representation of the generalized decoration-iteration transformation (DIT) 
for a more general decorating system composed of a finite number of decorating spins. Within the 
framework of this mapping method, the decorating system can be replaced with the effective pair 
interaction R between the two outer vertex spins a k \ and a k2 , as well as, the generally non- uniform 
magnetic fields H\ and H2 acting on those spins. 

loss of the validity of the algebraic transformation. For this most general case, the generalized 
decoration-iteration transformation can be formally written as 

Tr {Sfei} exp [-fiHk ({Ski}, Cfci, ^2)] = Aexp(f3Ra kl a k2 + /3Hia kl + f3H 2 a k2 ), (4.38) 

where the symbol Tr is h .\ denotes a trace over degrees of freedom of all decorating spins and 
the Hamiltonian 7-L k involves all the interaction terms depending on the decorating spins {S k i} 
and the two vertex Ising spins o k \ and <j k2 . It should be also mentioned that the aforelisted 
expressions f!4.36j) for the mapping parameters will remain valid whenever the Boltzmann's 
factor from the left-hand-side of the generalized decoration-iteration transformation (14.381) is 
used for obtaining the functions Vj (j — 1 — 4) instead of those given by the set of Eqs. ( 14.370 . 

4.3.2 Generalized star-triangle transformation 

Now, let us turn our attention to a generalization of the star-triangle transformation, which is 
applicable for the models in which each decorating system interacts with the three vertex Ising 
spins Ofci, cr k2 and a k3 . For simplicity, we will again consider first the decorating system being 
composed of a single decorating spin S k of the arbitrary magnitude and only then proceed to 
the subsequent generalization dealing with a more general decorating system being composed 
of a finite number of the decorating spins. The schematic representation of the generalized 
star-triangle transformation for the particular case of the decorating system consisting of a 
single decorating spin is displayed in Fig. [7] and this algebraic mapping transformation can be 

mathematically formulated as follows 

s 

ex P \fiS k (JiOfci + J 2 a k2 + J 3 a k2 )} 
s k =-s 
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G 



ki 




c STT, 




Figure 7: A diagrammatic representation of the generalized star-triangle transformation (STT) for the 
decorating system composed of a single decorating spin Sk of arbitrary magnitude. The terms Ji, J 2 
and J3 stand for the pair interactions between the decorating spin Sk and its three nearest-neighbour 
vertex spins a^i, &k2 and <7fc3> respectively. 

= A exp {PRia k2 (Jkz + PRi^kz^ki + (3R 3 ^ki^k2) • (4.39) 

The expression on the left-hand-side of the generalized star-triangle transformation ( 14 . 3 9 1) rep- 
resents the Boltzmann's factor of the four-spin cluster (star) from the left-hand-side of Fig. [7J 
which enters into the partition function of the Ising model taking into consideration three dif- 
ferent pair interactions J%, J 2 and J3 between the decorating spin Sk from a central site of the 
star and the three vertex spins <Jki, <Jk2 and <Jk?> from outer sites of the star. This expression is 
replaced via the generalized star-triangle transformation (I4.39P with the multiplicative factor A 
and the Boltzmann's factor of the three-spin cluster from the right-hand-side of Fig. which 
takes into account the pair interactions R\, R2 and R3 between the three outer vertex spins 
Cfei, 0fc2 and <7fc3 forming a triangle. The physical meaning of the star-triangle transformation 
( I4.39P thus lies in removing all the interaction terms associated with the central decorating 
spin Sk and substituting them through the effective pairwise interactions between the three 
outer vertex spins &ki, Ok2 and Ok?,- Of course, the generalized star-triangle transformation 
( I4.39P must hold independently of possible spin states of the three outer vertex spins and this 
self-consistency condition unambiguously determines the mapping parameters 
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which are expressed for the sake of brevity using the newly defined functions 
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Figure 8: A diagrammatic representation of the generalized star-triangle transformation (STT) for 
arbitrary decorating system composed of a finite number of the decorating spins. Within the framework 
of this mapping method, the decorating system can be replaced with the effective pair interactions 
Rl, R2 and R 3 between the three outer vertex Ising spins a k i, o k2 and cr k3 . 
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(4.41) 



Note furthermore that the functions Wj ( j = 1 — 4) label four different Boitzmann's weights, 
which can be obtained from the Boitzmann's factor on the left-hand- side of the generalized 
star-triangle transformation ( 14.391) by considering eight available spin configurations of the 
three vertex Ising spins (Tki, °"fc2 and 0-^3- 

Suppose now a more general decorating system, which comprises a finite number of the 
decorating spins rather than a single decorating spin, as it is schematically shown in Fig. El 
For this more general case, the generalized star-triangle transformation reads 



Tr {Ski} exp [-/3H k ({S ki }, a kXl a k2 , a k3 )} 

= A exp d3Rxa k2 a k 3 + P>Rio kz a kx + f3R 3 a kl a k2 ) , (4.42) 

where the symbol Tr^ Sk ^ denotes a trace over degrees of freedom of all decorating spins and 
the Hamiltonian 7-L k involves all the interaction terms depending on the decorating spins {S k i} 
and the three vertex Ising spins a k \, a k2 and a k3 . It should be mentioned that the set of 
expressions (14 .40 p . which is listed above for the relevant mapping parameters, remains valid if 
the Boitzmann's factor from the left-hand-side of the generalized star-triangle transformation 
(I4.42p is used for calculating the parameters Wj (j = 1 — 4) instead of those given by the set 
of Eqs. fOTj) . 
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4.3.3 On the validity of generalized transformations 

In this concluding part, let us briefly comment on a validity of the generalized algebraic trans- 
formations, which have been previously described in detail on two particular examples of the 
decoration-iteration and star-triangle mapping transformation. It becomes quite clear from a 
comparison of two preceding parts that the generalized decoration-iteration transformations 
(14.351) and (I4.38P hold in general even in a presence of the non-zero external magnetic field, 
while the generalized star-triangle transformations ( I4.39P and (14.421) hold true just in an absence 
of the external magnetic field. 

To clarify this fundamental difference, let us at first seek a general criterion determining 
a validity of the algebraic mapping transformations in a presence of the non-zero external 
magnetic field. If some general decorating system interacts with the m vertex Ising spins <jj, 
then, the total number of mapping parameters must be either greater than or at least equal 
to the total number of spin configurations available to the m vertex Ising spins, i.e., 2 m . It 
is worth mentioning that one mapping parameter can always be chosen in the form of the 
multiplicative factor, other m mapping parameters can be introduced by assuming the effect 
of inhomogeneous external magnetic field acting on the vertex spins and finally, m(m — l)/2 
mapping parameters can be ascribed to different pair interactions between the vertex spins. 
The generalized algebraic transformations in a presence of the external magnetic field should 
therefore obey the condition 



which is satisfied if and only if the decorating system interacts with either one or two outer 
vertex Ising spins (i.e. m < 2). The aforementioned criterion thus confirms a general validity 
of the decoration-iteration mapping transformation (m = 2) even in a presence of the non-zero 
external magnetic field, while the generalized star-triangle transformation (m = 3) apparently 
does not satisfy the demanded inequality (14.431) . 

On the other hand, the total number of independent configurations of the vertex Ising spins 
can be significantly lowered in an absence of the external magnetic field, because the relevant 
spin system has under this condition the spin inversion symmetry as it becomes invariant with 
respect to the reversal of all vertex Ising spins Oi — > —Oi (for 2 = 1,2,..., m). However, it should 




(4.43) 
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be emphasized that the total number of mapping parameters also reduces by canceling the m 
mapping parameters, which would represent the effect of inhomogeneous magnetic field acting 
on the vertex spins. Accordingly, the generalized algebraic transformations to be formulated in 
an absence of the external magnetic field should obey the inequality 



which holds true if and only if the decorating system interacts with either one, two or three 
outer vertex Ising spins (i.e. m < 3). This means that the generalized decoration- iteration and 
star-triangle transformations are valid under this constraint in accordance with the rigorous 
results reported on previously. Finally, it is worthy of notice that the inequalities (I4.43p and 
( 14.441) might serve in evidence that algebraic mapping transformations for any decorating system 
coupled to four or more outer vertex Ising spins cannot hold generally neither in a presence 
nor in an absence of the external magnetic field. This assertion of course holds true provided 
that only the pair spin-spin and single-spin interaction terms are taken into account when 
constructing the generalized mapping transformations. 

4.4 Generalized transformations II 

Another important progress in the development of generalized algebraic transformations has 
recently been achieved by Rojas, Valverde and de Souza [163] when taking into consideration 
a possibility of including higher-order multispin interactions in algebraic mapping transfor- 
mations with the aim to preserve their general validity. It should be noted here that under 
this circumstance it is always possible to find the algebraic mapping transformation of a quite 
general validity irrespective of the number of the vertex Ising spins involved in the interac- 
tions with the decorating system. Before constructing the most general star-polygon mapping 
transformation, the rigorous approach developed by Rojas et al. |163j will be adopted to treat 
a quite general decorating spin system interacting with either two or three vertex Ising spins 
within the framework of the generalized decoration-iteration and star-triangle transformation, 
respectively. 




(4.44) 
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4.4.1 Generalized decoration-iteration transformation 

Let us consider first arbitrary decorating spin system, which may be even of a quantum nature, 
coupled to the two outer vertex Ising spins a^i and <Jk2 as it is schematically shown in Fig. EJ 
The most general Hamiltonian of the spin system, which consists of all decorating spins and 
two vertex Ising spins, can be written as 

({Ski}, Cjfclj a k2) — ~ Jo'fi ({Ski}) — J\<Q ({Ski}) °k\ — Jq'} ({Ski}) °"fc2 

- 4 2 h{S kl })a k ia k2 . (4.45) 

Here, the relevant superscript marks a total number of outer vertex spins coupled with the 
decorating system and the symbol {S^} denotes a set of all decorating spins. The physical 
meaning of the interaction parameters Jj$ n ~ ( n i = 0, 1 for i = 1,2) is quite obvious; the 
former (latter) subscript m (n 2 ) determines whether or not the individual vertex spin <7ki 
(c"fc2) is indispensable part of the relevant interaction terms. Accordingly, the parameter Jq 2 q 
incorporates all the interaction terms depending only on the decorating spins {S^}. Among the 
most common examples of this type one could mention the interaction terms representing the 
effect of external magnetic field and zero-field splitting parameters on the decorating spins, as 
well as, bilinear, biquadratic and other higher-order interaction terms between the decorating 

(2) 

spins. Furthermore, the parameter J\ involves all the interaction terms, which depend on the 
particular vertex spin tiki and eventually on some of the decorating spins {Ski}- This parameter 
may for instance represent the influence of external magnetic field on the individual vertex spin 
(Tfci, as well as, the pair and higher-order interactions involving the vertex spin a^i and some of 
the decorating spins. Similar interaction terms are also included in the parameter Jq { except 
that the vertex Ising spin a k 2 now enters into the relevant interaction terms instead of a^x- 

(2) 

Finally, the parameter J\ { involves all the interaction terms depending on the product (Jki&k2 
between both vertex Ising spins and eventually on some of the decorating spins. In this respect, 
the pair interaction between the two outer vertex spins a k \ and <Jk2, as well as, the higher-order 
multispin interactions incorporating both these vertex spins and some of the decorating spins, 
are possible representatives of this type. 

Before presenting the most general form of the decoration-iteration transformation, it is 
unavoidable to perform a trace over degrees of freedom of all decorating spins in order to get 
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the effective Boltzmann's factor W k 2 \a k \, <J k2 ) depending merely on the two outer vertex Ising 
spins Ofci and o k2 



Wf Vfei, 0*2) = Tr {5fci} exp [-fiU^' ({S ki }, a ku a k2 )\ . (4.46) 

It should be remarked that the former step, which bears a close relation to capability of finding 
the relevant trace by undertaking analytical calculation (e.g. by exact analytical diagonaliza- 
tion), represents perhaps the most crucial limitation to applicability of the present method to 
any decorating system of a quantum nature. In the latter step, the effective Boltzmann's factor 
(14.461) can be replaced with a simpler equivalent expression depending just on the two outer 
vertex spins, which is provided by the generalized decoration-iteration mapping transformation 



,(2) 



Ki, cr k2 ) = exp (/3R£> + (3R\%a kl + (3R^{a k2 + (3R?ja kl a k2 ) . (4.47) 

It is of great practical importance that the self-consistency condition of the generalized decoration- 
iteration transformation (14.471) . which ensures a general validity of this mapping transformation 
irrespective of four possible spin states of the two vertex Ising spins involved therein, allows 
one to express all four yet undetermined mapping parameters through the unique formula 

W$, n2 = 2 2 (-+-)~ 2 £ £ al\all In [wf (a ku a k2 )} . (4.48) 

(2) 

The physical meaning of the relevant mapping parameters is as follows. The parameter R\ { 

stands for the effective pair interaction between the two outer vertex Ising spins, the parameters 

(2) (2) 

R[ and R { represent a generally non-uniform effective magnetic field acting on the vertex 

(2) 

spins (Tfci and a k2 , respectively, and the last parameter R q has merely a character of the 
multiplicative factor. With the help of the generalized decoration-iteration transformation 
given by Eqs. (I4.46P and (I4.47j) . the Hamiltonian (I4.45P of any decorating system coupled to 
two outer vertex Ising spins is effectively mapped onto a very simple Hamiltonian describing 
two mutually interacting Ising spins 

W Ki, a k2 ) = -R$ - i?gS<r fcl - R {2 }a k2 - R {2 }a kl a k2 . (4.49) 

Of course, the temperature-dependent effective interactions Rn} n2 °f corresponding model 
must obey the formula (I4.48P stemming from the self-consistency condition of the generalized 
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decoration-iteration transformation in order to validate the accurate mapping equivalence be- 
tween both Hamiltonians ( l4.45j) and (14.491) . Finally, it is also worthy to mention that the 
generalized form of the decoration-iteration transformation given by Eqs. (I4.46P and (I4.47P is 
fully consistent with the rigorous mapping transformation ( I4.38p . which has been previously 
proposed following the approach due to Fisher [162] . 

4.4.2 Generalized star-triangle transformation 

Even more interesting situation occurs in the search for the generalized star-triangle transfor- 
mation, which could be in principle applied to any decorating system coupled to the three outer 
vertex Ising spins as it is schematically illustrated in Fig. |HJ It has been argued previously that 
the rigorous approach developed by Fisher [162] does not establish the star-triangle transfor- 
mation of a quite general validity due to a lack of the mapping parameters, which represent 
all possible pair spin-spin and single-spin interactions. However, the single missing mapping 
parameter can be chosen so as to represent the effective triplet interaction between the three 
outer vertex spins when including this higher-order multispin interaction into a definition of 
the star-triangle transformation. 

The entire system, which is composed of all decorating spins and the three vertex Ising 
spins, can be described through the most general Hamiltonian of the following form 

({Ski}> <J k2, Cfe3) — — ^OAO ({Ski}) ~ •A^O ({Ski}) & kl ~ ^ofl,0 ({Ski}) °k2 

~ <^ofd,l ({Ski}) 0"fc3 — «^lfl,0 ({Ski}) Ok\Ok2 ~ «A)fl,l ({Ski}) &k2<?k3 

~ A%1 ({Ski}) 0"fc3CTfcl - ^1,1,1 ({Ski}) CTfcl0"fc20-fc3- (4.50) 

The physical meaning of the interaction parameters J$ n2 n3 (n, = 0, 1 for i — 1, 2, 3) is anal- 
ogous as before, namely, the notation suffix determines whether the individual vertex spin 
<Tfcj is present (n« = 1) or absent (rij = 0) in the interaction terms represented by the relevant 

(3) 

interaction parameter. For instance, the interaction parameter J| [ 1 involves all the interaction 
terms depending on the product Gki&k2&k3 of all three outer vertex Ising spins and eventually 
on some of the decorating spins {Ski}- 

After performing a trace over degrees of freedom of all decorating spins one obtains the 
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effective Boltzmann's factor 



Wjf\a kl ,a k2 ,(T k3 ) = Tr {5fci} exp [-(3U k A> ({S ki }, a kl , a k2 , a k3 )\ , (4.51) 

which depends solely on the three outer vertex Ising spins a k i, <*k2 an d cr k3 . The Boltzmann's 
factor (I4.5ip can be subsequently replaced with a simpler equivalent expression, which is sup- 
plied by the generalized star-triangle mapping transformation 

Wlf\a kl , a k2 , cr k3 ) = exp(/3i?JJ )0 + PR^ 0fi a kl + /3Rf} >0 a k2 + /^^i "^ 

+PRi] a kl a k2 + /3R^x!cr k 2(T k 3 + /3Rilicr m a kl + /3RU A(T k ia k2 a k3 ] . (4.52) 



,(3) 



The generalized star-triangle transformation (I4.52p satisfies the self-consistency condition, which 
demands its general validity regardless of eight available spin configurations of the three vertex 
Ising spins involved therein, if and only if the mapping parameters obey the unique formula 

PR$,n 2 , n3 = 2 2 (- + - + "3)-3 £ £ £ oj»jg In [wffai. * Mj a k3 )} . (4.53) 

Cfcl 0" fe2 cr fc3 

The physical meaning of the relevant mapping parameters is also quite evident. The mapping 

( 3 \ /q\ /q\ 

parameter Rq q is an appropriate multiplicative factor, while the parameters R\ q , -^0,1,0 an d 
(3) 

i?ooi represent a generally non- uniform effective magnetic field acting on the vertex spins a k i, 

(3) (3) (3) 

o k 2 and (Jfc3, respectively. Furthermore, the mapping parameters R\ { , -Rq 1 1 an d ^101 label 
the effective pair interactions between three exploitable couples of the vertex Ising spins and the 
parameter R\H denotes the effective triplet interaction among them. Using the generalized 
star-triangle transformation given by Eqs. (I4.5ip and (14.521) . the Hamiltonian (I4.50p of any 
decorating system coupled to the three vertex Ising spins is effectively mapped onto the most 
general Hamiltonian of three mutually interacting Ising spins 

-i>(3) / \ _ D (3) D (3) D (3) D (3) D (3) 

tt k [&kli CTk2,0' k3 ) — — -tCo q q — ri 10 Q<J k i — H 010 a k2 — -Ko.O,! ^ 3 — ftl t lfiO'klO'k2 

~ Ro,l,l tJ k20'k3 ~ R\ i fl,l°'k3< J kl — R\ i },l< J kl<3'k20'k3- (4.54) 

The rigorous mapping equivalence between both Hamiltonians (I4.50p and f)4.54p naturally de- 
mands that the temperature-dependent effective interactions Rn}n 2 n 3 °f the corresponding 
model are in concordance with the formula (I4.53P derived from a self-consistency condition of 
the generalized star-triangle transformation. 
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Figure 9: A diagrammatic representation of the generalized star-polygon transformation (SPT) for 
arbitrary decorating system composed of a finite number of the decorating spins. Within the framework 
of this mapping transformation, the decorating system coupled to arbitrary number of the outer 
vertex Ising spins can be replaced with the effective interactions, which are represented through edges 
connecting each couple of the vertex spins, as well as, any polygon whose vertices all coincide with 
lattice points of the vertex Ising spins. 

Finally, it should be nevertheless noted that an application of the generalized star-triangle 
transformation is of particular research interest just if it establishes a precise mapping equiva- 
lence with some simpler exactly solvable model. It is therefore valuable to mention few rigor- 
ously solved Ising models of this type. The spin-1/2 Ising model with only the triplet interaction 
is exactly tractable on the triangular lattice (the so-called Baxter- Wu model) [164H168] . the 
union jack lattice |169H171j and the diced lattice |172[I173] . Among the more general exactly 
solved Ising models extended by the triplet interaction one could also mention the spin-1/2 Ising 
model on the diced |174j and union jack |175] lattices with the pair and triplet interactions, 
the spin-1/2 Ising model on the union jack lattice with the triplet interaction and the external 
magnetic field |176H178j . as well as, the spin-1/2 Ising model on the kagome lattice |179j and 
Cayley tree [180] including the external magnetic field, the pair and triplet interactions. 



4.4.3 Generalized star-polygon transformation 

Last but not least, let us examine the most general form of the star-polygon transformation, 
which would enable to treat in principle any decorating system coupled to arbitrary number 
of the vertex Ising spins as it is schematically shown in Fig. |HJ It should be mentioned that 
the generalized star-polygon transformation involves the generalized decoration-iteration and 
star-triangle transformations as the special but surely the most notable cases. The most general 
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Hamiltonian for the entire spin system, which consists of all decorating spins and the m vertex 
Ising spins, can be written in the form 

m 

4 m) ({Su}, = -£ J$ } ({S ki }) II *3- ( 4 -55) 

{rij} 3=1 

Here, the symbol {rij} = {ni, n%, . . . , n m } denotes a set of all two-valued variables, each of 
them determining whether the individual vertex spin o^j is present (rij = 1) or absent (rij = 0) 
in the interaction terms represented by the interaction parameter Jj^y- Next, the sets {Ski} = 
{S k i, S k2 , S kn } and {o k j} = Wki, &k2, ■ ■ ■ , o'km} span all the decorating spins and vertex 
Ising spins, respectively, and the relevant summation £ = £ £ ' ' ' £ ^ s carr i e d out 

{rij} m=0,lri2=0,l n m =0,l 

over the full set of two- valued variables {rij}. 

The generalized star-polygon transformation can be introduced in three crucial steps. In 
the first step, it is necessary to perform a trace over degrees of freedom of all decorating spins 
in order to get the effective Boltzmann's factor 



W, 



(m). 



k '(Wkj}) = Tr {5fci} exp [-/3?4 m) ({S ki }, {a kj })\ , 

which is solely dependent just on the m vertex Ising spins from the set {<Jkj}- In the second step, 
the Boltzmann's factor (I4.56P can be substituted by a simpler equivalent expression provided 
by the generalized star-polygon mapping transformation 



(4.56) 



Wi m) (Ki})=ex P 



e ( n - 



.{%} 



3=1 



(4.57) 



In the third step, it is easy to prove that the generalized star-polygon transformation (I4.57P is 
valid independently of all spin configurations available to the m vertex Ising spins if and only 
if the relevant mapping parameters fulfill the unique formula 



i m 



In 



wt\Wk 3 }) 



(4.58) 



Bearing all this in mind, the generalized star-polygon transformation given by Eqs. (I4.56P 
and ( 14.57P establishes a rigorous mapping equivalence between the Hamiltonian ( 14. 55ft of any 
decorating system coupled to the set {dkj} of the vertex Ising spins and respectively, the most 
general Hamiltonian describing of the m mutually interacting vertex Ising spins including all 
higher-order multispin interactions 



H k 



{Wkj}) 



{n 3 } j=l 



(4.59) 
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The rigorous mapping equivalence between both Hamiltonians (14.551) and (14.591) of course holds 
just if the temperature-dependent effective interactions R^ y of the corresponding model meet 
the unique formula (I4.58p . which has been obtained from the self-consistency condition of the 
generalized star-polygon transformation. 

To conclude, the problem of finding a rigorous solution for any lattice-statistical model in 
which the decorating system is coupled to arbitrary number of the vertex Ising spins turns out 
to be equivalent with the problem of solving the corresponding Ising model with temperature- 
dependent higher-order multispin interactions between those vertex Ising spins. Another in- 
teresting observation is that the highest order of effective multispin interactions is in general 
equal to the total number of the vertex Ising spins to which the relevant decorating system is 
coupled. In the consequence of that, the generalized star-square transformation as applied to 
any decorating system interacting with four outer vertex spins will for instance afford a pre- 
cise mapping relationship with the respective Ising model taking into consideration multispin 
interactions up to the fourth order. From this point of view, the lattice-statistical models in 
which the decorating system is coupled to four or more outer vertex spins are just barely fully 
exactly solvable due to a lack of exactly solved Ising models accounting for higher-order mul- 
tispin interactions. There are however few notable exceptions among which one could mention 
the spin-1/2 Ising model with the pair interaction on two square lattices coupled together by 
means of the quartic interaction, which is equivalent with a general eight-vertex model as first 
evidenced by Wu |181j . Kadanoff and Wegner |182j . Albeit the general eight- vertex model is 
not exactly tractable, two valuable cases of the eight-vertex model has been rigorously solved 
under the special constraints to its Boltzmann's weights. The symmetric eight-vertex model 
satisfying the zero-field condition has exactly been solved by Baxter |183[I184] and the free- 
fermion eight-vertex model obeying the free-fermion condition has rigorously been solved by 
Fan and Wu |185[I186] . Among other matters, the utility of the generalized star-square trans- 
formation to provide exact results from a precise mapping equivalence with either the zero-field 
or free-fermion eight- vertex model will be illustrated on two particular exactly tractable models 
in the following two sections. 
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5 Exactly solved Ising models 

In this section, let us demonstrate the undeniable capability of generalized algebraic transfor- 
mations to enrich the realm of exactly solved Ising models. For this purpose, four scientific 
papers of the present author and collaborators are listed in the appendices A1-A4, which 
constitute the supplementary part of this book not included therein. These supplements are 
concerned with the Ising models, which have been exactly solved by taking advantage of the 
generalized decoration-iteration, star-triangle and star-square transformation, respectively. It 
is worth mentioning that the exactly solved Ising models, which have been earlier treated with 
the help of algebraic mapping transformations, are reviewed in the introductory parts of both 
papers from the appendices Al and A2. In what follows, let us briefly summarize the most 
important scientific achievements of the papers from the appendices A1-A4, whereas the reader 
interested in more details is referred to the relevant articles for more details. 

5.1 Decoration- iteration transformation 

In the first paper from the series on the exactly solved Ising models, which is presented as the 
Appendix Al, the generalized decoration- iteration transformation is employed to obtain the 
exact solution for the mixed-spin Ising model on a decorated square lattice with two different 
kinds of decorating spins placed on its horizontal and vertical bonds. At first sight, this exactly 
solvable model might seem as a rather trivial extension of the mixed-spin Ising model on a 
decorated square lattice with the same kind of decorating spins on all horizontal as well as 
vertical bonds, which has been proposed and exactly solved by Jascur |187j and Dakhama [188] 
almost a decade ago. The most essential difference between both models lies in the fact that 
the former model with two different types of decorating spins is effectively mapped onto the 
spin-1/2 Ising model on the anisotropic square (rectangular) lattice, while the latter model 
with just one kind of decorating spins is mapped onto the spin-1/2 Ising model on the isotropic 
square lattice. 

It should be nevertheless mentioned that this more or less irrelevant difference has a very 
significant impact on the diversity of available spontaneous long-range orderings. In addition 
to the classical ferromagnetic or ferrimagnetic spin arrangement, the mixed-spin Ising model 
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with the integer and half-odd-integer decorating spins placed on horizontal and vertical bonds 
of a square lattice may also exhibit a peculiar spontaneous long-range order of the quasi- ID 
character. This unusual spontaneous order appears as a result of the sufficiently strong (but 
not too strong) uniaxial single-ion anisotropy, which forces all integer-valued decorating spins 
towards their non-magnetic spin state. Owing to this fact, the mixed-spin Ising model on 
a decorated square lattice effectively splits into a set of quasi- independent mixed-spin Ising 
chains, which comprise the alternating vertex Ising spins and half-odd-integer decorating spins. 
However, the most important finding to emerge from this study is that the investigated model 
surprisingly remains spontaneously long-range ordered in some restricted range of the uniaxial 
single-ion anisotropics despite the obvious effective reduction of its spatial dimensionality. As 
a matter of fact, the critical temperature tends to zero just at some more negative uniaxial 
single-ion anisotropy as one would intuitively expect (see Fig. 3 in the Appendix Al) and 
moreover, the spontaneous sublattice magnetization of the integer-valued decorating spins is 
raised from zero just on behalf of selective thermal fluctuations (see Fig. 5(d) in the Appendix 
Al). These observations have obvious relevance to a deeper understanding of critical behaviour 
of quasi-lD spin systems prone to a spontaneous ordering below some critical temperature, 
which need not necessarily arise from much weaker interactions establishing spontaneous order 
of 3D character but may represent inherent feature of a relevant quasi-lD system instead. Of 
course, this outstanding feature cannot be observed in the mixed-spin Ising model with just 
one kind of the decorating spins, which has been explored by Jascur |187] and Dakhama [188 J 
in the earlier publications. 

5.2 Star-triangle transformation 

The mixed spin- 1/2 and spin-S* (S > 1/2) Ising model on a bathroom-tile lattice, which is 
the main objective of the second paper presented in the Appendix A2, has been rigorously 
treated by adopting the generalized star-triangle and triangle-star transformations. The most 
important focus of this work is to reveal the main difference between the magnetic behaviour 
of two different but relative mixed-spin Ising models on the regular honeycomb and semi- 
regular bathroom-tile lattice. It is worthwhile to remark that the honeycomb and bathroom-tile 
lattices are two planar lattices with the same coordination number three but a rather different 
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lattice geometry. The regular honeycomb lattice constitute congruent polygons of only one type 
(hexagons) and accordingly, all vertices as well as all edges of this 2D lattice are equivalent. 
Contrary to this, the semi-regular bathroom-tile lattice has equivalent only all vertices (but not 
edges), because this Archimedean lattice consists of two different regular polygons - squares 
and octagons [161] . 

It has been already mentioned before that a comparison between critical points of the 
spin-1/2 Ising model on the regular square and semi-regular kagome lattice (with the same 
coordination number four) has implied a slightly higher critical temperature of the former reg- 
ular lattice [16U] . The critical temperature apparently reflects a robustness of the spontaneous 
long-range order against thermal fluctuations and it might be therefore quite interesting to 
ascertain whether the planar lattices with the same coordination number but higher symmetry 
always have higher critical temperature than their less symmetric counterparts. Namely, it 
is tempting to conjecture that this behaviour is universal and holds true for any Ising model 
irrespective of the coordination number, spin magnitude, etc. The paper supplied in the Ap- 
pendix A2 provides convincing evidence that the same conclusion is reached when comparing 
the respective critical points of the spin-1/2 Ising model on the honeycomb and bathroom-tile 
lattice, since the critical temperature of the honeycomb lattice lies just slightly above that of 
the bathroom-tile lattice. Moreover, the more general mixed spin-1/2 and spin-S 1 Ising model 
on the honeycomb and bathroom-tile lattices also behaves completely in agreement with this 
assertion, which is evidently not affected neither by the higher spin magnitude S > 1, nor by 
the uniaxial single- ion anisotropy (see for instance Fig. 2 in the Appendix A2). 

Finally, it should be stressed that another valuable by-product of calculations presented 
in the paper from the Appendix A2 is the exact solution of the spin-1/2 Ising model on the 
Shastry-Sutherland (orthogonal-dimer) lattice, which has been published in the subsequent 
separate work in full details [189] . Even although the spin-1/2 Ising model on any planar lattice 
is in principle exactly tractable problem within the framework of Pfaffian technique [92H94"] . 
this rigorous method usually turns into a rather cumbersome and tedious calculation when 
applying it to some more complex planar lattice of lower symmetry. At the present state of 
knowledge, it seems quite striking that the critical temperature of the spin-1/2 Ising model 
has not been exactly evaluated yet for all semi-regular Archimedean lattices, which are the 
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only plane tilings having all vertices equivalent |161j . To the best of our knowledge, the Monte 
Carlo simulations has been implemented as the only numerical method in order to locate the 
respective critical points of the spin- 1/2 Ising model on all twelve Archimedean lattices [190J 
and thus, this analytically unresolved problem represents another challenging issue to deal with 
in the future by making use of generalized algebraic mapping transformations. 

5.3 Star- square transformation 

Another two papers provided in the Appendices A3 and A4 deal with the mixed spin-1/2 
and spin-S" Ising model on the union jack (centered square) lattice, which can be rigorously 
mapped onto the corresponding eight-vertex model with the help of generalized star-square 
transformation. It is noteworthy that the model under consideration represents an interesting 
example of geometrically frustrated spin system, which may exhibit on account of a competition 
between the nearest-neighbour and next-nearest-neighbour interaction a rather rich critical 
behaviour including reentrant phase transitions, non-universal critical behaviour, as well as, 
discontinuous first-order transitions. 

First, let us make few remarks on a critical behaviour of several particular limiting cases. 
The special case of the spin-1/2 Ising model on the union jack lattice is fully exactly solvable 
due to a precise mapping equivalence with the free-fermion eight-vertex model as convincingly 
evidenced by Wu and Lin [191]. In addition, this particular case simultaneously represents the 
first exactly solved Ising model displaying reentrant phase transitions in a relatively narrow 
region of the parameter space |192j . Note furthermore that this special case is otherwise less 
interesting, since its critical behaviour belongs to the standard Ising universality class and be- 
sides, this model does not even exhibit first-order phase transitions. The much more interesting 
situation emerges in the particular cases with the higher spins S > 1, which are not completely 
exactly soluble on behalf of a rigorous mapping equivalence with the more general eight-vertex 
model. Exact results for those particular cases are therefore restricted only to a certain manifold 
of the parameter space, which ensures a precise mapping correspondence either with the zero- 
field |183|J184j or free-fermion [185,186j eight- vertex model. In the rest of the parameter space, 
the rather meaningful and reliable results can be obtained from the free-fermion approximation 
simply ignoring a non- validity of the free-fermion condition |185|,I186] . 
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It has been proved by Lipowski and Horiguchi [193] that the mixed spin- 1/2 and spin- 
1 Ising model on the union jack lattice exhibits a striking non- universal critical behaviour 
in addition to reentrant phase transitions. The non-universal criticality is closely connected 
with an appearance of special multicritical points usually referred to as bicritical points, at 
which two second-order and one first-order phase transition lines meet together. Moreover, 
Lipowski and Horiguchi |193] have furnished a rigorous evidence that the critical exponents 
vary continuously along the line of bicritical points in accordance with the weak universality 
hypothesis proposed by Suzuki [194] . It should be emphasized that bicritical points as well as 
first-order phase transitions have been calculated from a precise mapping relationship with the 
zero-field eight-vertex model and thus, they cannot be artefacts of any approximation because 
of their exactness. The essentially same critical behaviour has been also detected in the mixed 
spin-1/2 and spin-3/2 Ising model on the union jack lattice, which has been examined by the 
present author in the publication included in the Appendix A3. However, the most fundamental 
difference between both aforementioned mixed-spin Ising models has been found in the allowable 
values of the critical exponents, which continuously vary along the line of bicritical points. In 
fact, the availability range of the critical exponents for the mixed-spin Ising model on the union 
jack lattice with S = 1 has turned out to be roughly twice as large as that calculated for the 
analogous mixed-spin Ising model with S = 3/2. 

To clarify this intriguing issue, another two particular cases of the mixed spin-1/2 and 
spin-S 1 Ising model on the union jack lattice with S = 2 and S = 5/2 have been explored 
in the comprehensive work listed in the Appendix A4. This theoretical study has confirmed 
correctness of conclusion on two different ranges of availability of the critical exponents, which 
are governed by the specific condition of whether the higher spin S is integer or half-odd-integer 
one. 
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6 Exactly solved Ising-Heisenberg models 

In this section, let us demonstrate a versatility of the generalized algebraic transformations in 
providing exact results for the lattice-statistical models in which the decorating system is chosen 
in the form of a rather small cluster of the quantum Heisenberg spins. The hybrid classical- 
quantum spin systems composed of the semi-classical Ising spins and the quantum Heisenberg 
spins will be consequently described within the framework of Ising-Heisenberg models, which 
are in principle exactly tractable by the use of rigorous approach developed on the grounds 
of generalized decoration-iteration, star-triangle and star-square transformations, respectively. 
Before proceeding to a short survey of the most notable findings, which are reported in detail 
in the papers from the Appendices B1-B4 dealing with three exactly solved spin- 1/2 Ising- 
Heisenberg models, it should be pointed out that new exciting findings might be expected on 
account of a presence local quantum fluctuations. From this perspective, it is of particular 
research interest to explore especially 2D Ising-Heisenberg models, which might simultaneously 
exhibit a rather unusual combination of both spontaneous long-range order as well as obvious 
quantum features. 

6.1 Decoration- iteration transformation 

The spin-1/2 Ising-Heisenberg model on doubly decorated planar lattices, which is the main 
subject of the article from the Appendix Bl, has been accurately solved by means of the 
generalized decoration-iteration transformation. The most rewarding aspect of this rigorous 
study surely represents a discovery of unconventional antiferromagnetic spontaneous long-range 
order, which occurs in spite of a pure ferromagnetic character of all considered interactions. 
It should be mentioned that this striking and rather unexpected quantum antiferromagnetic 
phase (QAP) appears as a result of the mutual competition between two different ferromagnetic 
interactions; the former is the Ising interaction J\ of the easy-axis type and the latter one is 
the XXZ Heisenberg interaction J (A) of the easy-plane type. Indeed, it has been proved that 
the classical ferromagnetic order in which all Ising as well as Heisenberg spins are aligned 
parallel one to each other represents the lowest-energy state (ground state) just if the exchange 
anisotropy A in the XXZ Heisenberg interaction is less than the critical value A c = y / T+ 2 Ji/J, 
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while the unconventional QAP becomes the ground state in the reverse case A > A c > 1. 

Let us take a closer look at the QAP since its mere appearance is a very surprising manifes- 
tation, which cannot be evidently of a classical nature but of a quantum origin. For simplicity, 
the unusual spin order emerging in the QAP will be described for the particular case of the 
spin-1/2 Ising-Heisenberg model on doubly decorated square lattice (see Fig. 1 in the Appendix 
Bl), where it can be expressed through the eigenvector written in the following compact and 
suggestive form 



.lY+j 



1 



iqap) = n i > n 2' - 2> + sin0 (fU 



1 1 



2 2l>T> 



Here, the mixing angle is given by the formula <j> = | arctan f^J , the former product runs 
over all nodal sites of a square lattice occupied by the vertex Ising spins, the indices i and j 
label the relevant row and column of a square lattice, respectively, and the latter product is 
carried out over all bonds of a square lattice containing the Heisenberg spin pairs (dimers). The 
eigenvector (16 .ip indicates that the perfect Neel order is captured on the nodal sites of a square 
lattice by the vertex Ising spins, while the respective behaviour of the Heisenberg dimers is 
governed by the quantum entanglement of two antiferromagnetic microstates available to each 
couple of the Heisenberg spins. It should be noticed, however, that the mutual orientation 
of the nearest-neighbour Ising and Heisenberg spins is neither perfectly parallel nor perfectly 
antiparallel on behalf of a zero-point motion of the Heisenberg spins, which is triggered by 
the local quantum fluctuations. According to Eq. (16.11) . the parallel (antiparallel) alignment 
of the nearest-neighbour Ising and Heisenberg spins always represents the more (less) probable 
microstate occurring with the probability amplitude cos</> (sin0), which gradually diminishes 
(enhances) upon strengthening the ratio JA/J^ It is noteworthy that the aforedescribed nature 
of the QAP is consistent with an independent analysis of spin dynamics, which has been carried 
out by exploring time-dependent autocorrelation and correlation functions [195J. 

As far as the finite-temperature behaviour is concerned, the model under investigation ex- 
hibits quite interesting global phase diagram in which two second-order phase transition lines 
separating the classical ferromagnetic phase and the unconventional QAP from the param- 
agnetic region merge together at a quantum critical point (see Fig. 5 in the Appendix Bl). 
It should be also noted here that both critical lines are from the standard Ising universality 
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class as convincingly evidenced by a detailed analysis of the relevant order parameters, which 
is the spontaneous magnetization in the classical ferromagnetic phase and the staggered mag- 
netization in the unusual QAP. The same conclusion can be found from the detailed analysis 
of specific heat, which displays at a critical point of the order- disorder transition a common 
logarithmic singularity regardless of whether the classical ferromagnetic or the unconventional 
QAP constitutes the ground state (see Fig. 5 in the Appendix Bl). In future, it could be quite 
interesting to ascertain whether an existence of the unconventional QAP represents a generic 
feature of wider family of quantum spin models with competing easy-axis and easy-plane fer- 
romagnetic interactions or it is just a specific feature of the hybrid Ising-Heisenberg models. 
The exactly solved Ising-Heisenberg models could provide a useful starting basis for performing 
perturbative calculations, which could bring insight into this unresolved problem. 

6.2 Star-triangle transformation 

The generalized star-triangle transformation has been also adapted in order to obtain the ex- 
act solution for the spin- 1/2 Ising-Heisenberg model on the triangulated kagome (triangles- in- 
triangles) lattice, which is supplemented in the paper listed in the Appendix B2. It is worthy 
to notice that the geometric structure of the model under investigation has been motivated 
by the magnetic lattice of a series of three isostructural polymeric coordination compounds 
Cu9X 2 (cpa)6.nH 2 (X = F, CI, Br and cpa = carboxypentonic acid) |196H198| . The magnetic 
structure of this family of transition-metal complexes displays a rather rare and curious archi- 
tecture, which is constituted by smaller triangular-shaped spin clusters embedded in greater 
triangular unit cells forming the kagome pattern |196H198] . In addition, the strong antifer- 
romagnetic exchange interaction between the decorating spins from smaller triangular-shaped 
entities has been evidenced by several experimental studies [199H2U2] and thus, this family of 
magnetic compounds belongs to a prominent class of highly frustrated magnetic materials. In 
accordance with this statement, all three isostructural magnetic compounds from this family 
does not order down to 1.3K despite the impressive exchange interactions of few hundreds of 
Kelvin [199T I2U2] and this finding might be regarded as experimental indication of a disordered 
spin-liquid ground state. 

With this background, the primary goal of the article presented in the Appendix B2 has 

51 



6.3 Star-square transformation 



6 ISING-HEISENBERG MODELS 



been focused on the role of local quantum fluctuations in determining the nature of disordered 
spin-liquid state. For this purpose, we have precisely examined both spin- 1/2 Ising-Heisenberg 
model on the triangulated kagome lattice as well as its analogous semi-classical Ising model 
with the aim to shed light on differences in their respective behaviour, which might be very 
pronounced due to a presence (absence) of local quantum fluctuations in the former (latter) 
model. Our exact calculations have furnished a rigorous proof that the residual entropy of 
the spin-liquid phase inherent to the Ising-Heisenberg model with any but non-zero exchange 
anisotropy amounts roughly a half of that calculated for the analogous classical Ising model. 
Owing to this fact, the local quantum fluctuations basically diminish a macroscopic degeneracy 
of the disordered spin-liquid state even though the residual entropy is not completely removed 
by the quantum ' or der-from- disorder' effect unless the antisymmetric Dzyaloshinskii-Moriya 
interaction is considered [203J. 

Besides a comprehensive analysis of the disordered spin-liquid state, the critical behaviour 
and basic thermodynamic quantities of the spin-1/2 Ising-Heisenberg model on the triangu- 
lated kagome lattice have been also investigated in dependence on a strength of the geometric 
frustration. It has been evidenced that the stronger the exchange anisotropy in the XXZ Heisen- 
berg interaction is, the weaker antiferromagnetic interaction is needed in order to destroy the 
spontaneous long-range order. In addition, the critical temperature exhibits a striking non- 
monotonous dependence on a strength of the geometric frustration whenever the easy-plane 
XXZ Heisenberg interaction competes with the easy-axis Ising interaction. It is worthy of 
notice, moreover, that our exact calculation has also enabled us to conjecture the notable 
temperature dependence of the specific heat with two well separated round maxima, which 
should be experimentally observed in the three isostructural polymeric coordination compounds 
Cu9X 2 (cpa)6.nH 2 0. Finally, it is worthwhile to remark that several exact results presented in 
the Appendix B2 have been confirmed by a subsequent independent calculation performed by 
Yao et al. [201] . 

6.3 Star- square transformation 

The last two papers provided in the Appendices B3 and B4 deal with the spin-1/2 Ising- 
Heisenberg model with the pair XYZ Heisenberg interaction and the quartic Ising interaction. 
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It should be pointed out that the model under investigation falls into a prominent class of fully 
exactly solved models, which are rigorously tractable due to a precise mapping equivalence 
with the zero-field eight-vertex model established through the generalized star-square transfor- 
mation. The relevance of the exactly solved Ising-Heisenberg model with the pair Heisenberg 
interaction and the quartic Ising interaction should be viewed in more respects. The most 
attractive issue of this exactly solved model certainly represents a violation of the strong uni- 
versality hypothesis, which could be intuitively expected already from a rigorous mapping 
relation with the zero-field eight-vertex model. However, the specific details of how critical 
exponents depend on relevant interaction parameters remains obscure until the comprehensive 
analysis is made to resolve availability range of the critical exponents. 

First, it is worthwhile to remark that the model with the antiferromagnetic pair interaction 
surprisingly exhibits less evident changes of both critical temperatures as well as critical expo- 
nents than the model with the ferromagnetic pair interaction and henceforth, our attention will 
be mainly concentrated on the latter model. The critical temperature of the spin-1/2 Ising- 
Heisenberg model with the ferromagnetic XXZ Heisenberg interaction displays a remarkable 
dependence on the exchange anisotropy with two marked wings of critical lines, which merge 
together at a quantum critical point between two different spontaneously long-range ordered 
phases (see Fig. 4 in the Appendix B3 and Figs. 2-4 in the Appendix B4). It should be em- 
phasized that this peculiar quantum critical point is accompanied with a singular behaviour of 
the critical exponents and consequently, it bears evidence of a phase transition of the infinite 
order (see Fig. 6 in the Appendix B3 and Figs. 5-7 in the Appendix B4). If the XXZ ferro- 
magnetic interaction is assumed, then, reentrant phase transitions may be also observed in a 
close vicinity of the quantum critical point due to different degeneracies of both spontaneously 
long-range ordered phases. However, the observed reentrance is gradually suppressed by con- 
sidering the less symmetric XYZ ferromagnetic pair interaction as convincingly evidenced in 
Figs. 5-7 depicted in the Appendix B4. It has been proved that the disappearance of reentrant 
phase transitions really occurs because of the more general XYZ exchange anisotropy, which 
generally lifts a macroscopic degeneracy of the one from two spontaneously long-range ordered 
phases. 

To summarize, the spin-1/2 Ising-Heisenberg model with the pair XYZ Heisenberg and 
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quartic Ising interactions represents a rare example of the exactly solved classical-quantum 
model with a striking weak-universal critical behaviour. It has been demonstrated that the 
critical exponents may continuously vary over the entire range of allowable values when the 
relevant interaction parameters are appropriately chosen. Several interesting extensions and 
generalizations of the present version of the Ising-Heisenberg model come also into question. For 
instance, it is possible to extend the present model by including higher-order triplet, quintuplet 
and sextuplet interactions between the Ising and Heisenberg spins or to study thermodynamic 
behaviour of this exactly solved model in more detail. This might serve as motivation for future 
work. 
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7 Conclusions and future outlooks 

The present book concerns with the theory of generalized algebraic transformations, which are 
capable of producing new interesting exactly solvable models after performing relatively modest 
calculations. This implies a considerable relevance of algebraic mapping transformations, since 
the exactly soluble models of interacting many-body systems are currently considered as an 
inspiring research field in its own right [IHT3]. As a matter of fact, the search for exactly solv- 
able models is usually recognized as the highest intellectual challenge for theoretical physicists 
due to a rather intricate and sophisticated mathematics, which is usually involved in any at- 
tempt to treat even relatively simple lattice-statistical models rigorously [205|, I206] . From this 
viewpoint, the rigorous method based on generalized algebraic transformations avoids those 
formidable mathematical difficulties by establishing a precise mapping equivalence with some 
simpler lattice-statistical models whose exact solutions are already known. Notwithstanding 
that this work has been aimed at promoting the implementation of generalized algebraic trans- 
formations, which can be adapted to discover completely novel and yet unexplored phenomena 
without a direct resemblance with any observable phenomenon of the corresponding exactly 
solved model. 

The importance of generalized algebraic transformations has been convincingly evidenced 
in two wide families of the exactly solvable models. The series of exactly solved Ising models 
brought a deeper insight into diverse aspects closely connected especially with phase transitions 
and critical phenomena. More specifically, the exactly solved Ising models from the Appendices 
A1-A4 shed light on a critical behaviour of spontaneously long-range ordered quasi-lD systems, 
reentrant phase transitions, non-universal critical behaviour, magneto-structural correlations 
and so on. On the other hand, the exactly solved Ising- Heisenberg models from the Appendices 
B1-B4 have served in evidence that the generalized algebraic transformations can be even 
adapted to treat accurately the hybrid classical-quantum models as well. In this class of the 
exactly solved models, our attention has been primarily focused on whether spontaneous long- 
range order might be accompanied with obvious macroscopic features of a quantum origin. 
Among the most remarkable findings, which have been reported on the exactly solved Ising- 
Heisenberg models, one could mention the unconventional antiferromagnetic spontaneous long- 
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range order arising from a competition between two different ferromagnetic interactions, the 
rigorous analysis of quantum critical points and weak-universal critical behaviour, as well as, the 
partial lifting of a macroscopic degeneracy in disordered spin liquid states owing to a presence 
of local quantum fluctuations. 

Before concluding, it is worthwhile to remark that the rigorous technique based on the 
grounds of generalized algebraic transformations offers a variety of other opportunities to deal 
with in the future and the presented exactly solved models hopefully illustrate the most fun- 
damental aspects of this powerful exact mapping method. It has been already argued in a 
breakthrough article published more than a half century ago by Fisher that any statistical- 
mechanical system might be selected as the decorating system without disturbing the validity 
of generalized algebraic transformations [162J. Obviously, this fact implies that the generalized 
algebraic transformations might be employed for diverse lattice-statistical models including 
the ones of hybrid classical-quantum nature. Apart from the exactly solved Ising-Heisenberg 
models, the algebraic mapping transformations have been recently engaged also for another 
intriguing class of hybrid classical-quantum models describing the mutually interacting system 
composed of the localized Ising spins and delocalized electrons [207H210] . It is the author's 
hope that many challenging and so far unresolved issues in the area of exactly solvable models 
will be tackled in the future within the framework of this rigorous mapping method. 
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